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EXECUTIVE  PROGRAM  SUMMARY 


This  report  describes  the  x-ray  laser  research  conducted  by  the  Radiation 
Hydrodynamics  Branch  under  SDIO  sponsorship  during  FY  1992.  Some  of  the 
experimental  work  was  performed  at  Sandia  National  Laboratories  as  part  of  a 
collaboration  with  NRL.  This  document  consists  of  four  sections  whose  authors  are 
separately  identified,  each  detailing  progress  in  a  different  aspect  of  the  overall  effort 

Last  year’s  report  described  the  first,  and  still  the  only,  achievement  of  a  population 
inversion  in  the  x-ray  region  of  the  spectrum  by  resonant  photopumping.  The  experiment 
consisted  of  two  basic  components:  an  intensely  radiating  pure  sodium  Z  pinch  whose 
He  a  transition  is  resonant  at  11  A  with  the  1-4  He  7  line  of  neon,  and  a  neon  gas  cell 
which  is  photoionized  by  the  approximately  10  TW  of  total  x-rays  coming  from  the  pinch. 
Since  that  noteworthy  accomplishment,  further  work  has  concentrated  on  enhancing  the 
photoionization  of  the  neon  and  detecting  and  measuring  the  possible  laser  lines. 

To  photoionize  neon  to  the  He-like  stage,  it  is  critical  that  the  pumping  pinch  produce 
sufficient  photons  in  the  energy  regions  near  the  photoionization  edges  of  the  Li,  Be,  and 
B-like  stages  (about  250  -  450  eV  for  Ne).  The  first  section  describes  work  in  which 
the  soft  power  output  spectrum  of  the  pure  sodium  pinch  was  measured  as  a  function 
of  energy.  These  initial  measurements  showed  that  the  photon  flux  in  the  critical  250  - 
450  eV  region  was  considerably  less  than  that  which  would  obtain  if  the  pinch  radiated  a 
Planckian  spectrum  with  a  temperature  appropriate  to  the  total  measured  power  output. 
As  a  means  of  increasing  this  critical  photon  power,  potassium  was  introduced  into  the 
(previously)  pure  sodium  pinch,  since  it  is  a  more  efficient  soft  x-ray  radiator.  The  results 
are  easily  summarized:  it  works.  Replacing  25%  of  the  sodium  with  potassium  has  no 
appreciable  effect  on  the  power  of  the  optically  thick  pump  line,  but  raises  the  subkilovolt 
power  of  the  pinch  a  remarkable  factor  of  3.  Furthermore,  the  comparative  spectra  of  neon 
clearly  demonstrate  enhanced  ionization  when  the  sodium-potassium  load  is  employed. 

The  figure  at  the  end  of  this  summary  shows  three  of  the  spectra  in  the  region  of  the 
possible  neon  laser  lines  which  were  obtained  at  Sandia  with  the  new  sodium-potassium 
load.  The  two  lines  near  217  A  may  represent  the  first  detection  of  a  He-like  4-3  line  but 


this  identification  is  still  uncertain  at  this  time.  Note  also  the  unexplained  increase  of  a 
factor  of  5  of  the  5f-3d  Li-like  line  when  the  fill  pressure  is  increased  by  only  50%.  The 
5f-3d  line  in  Li-like  A1  has  been  observed  to  lase  in  recombination  by  no  less  than  8  groups 
worldwide.  This  very  interesting  observation  in  our  photopumped  system  is  now  under 
further  experimental  and  theoretical  investigation  for  the  potential  of  a  recombination 
laser. 

Another  important  advance  in  photopumped  x-ray  laser  research  was  accomplished  at 
the  Laboratory  of  Plasma  Studies,  Cornell  University,  in  work  which  had  been  supported 
in  part  on  a  previous  SDIO  contract  and  had  involved  some  collaboration  with  NRL.  At 
Cornell,  resonant  photopumping  has  been  demonstrated  via  fluorescence  observations  in 
the  A1  -  Mg  system,  on  a  pulsed-power  device  producing  only  400  kiloamperes  peak 
current  This  work  has  been  accepted  for  publication  by  Physical  Review  A  and  is 
scheduled  for  the  March,  1993  issue.  The  authors  are  Qi,  Hammer,  Kalantar,  and  Mittal. 

Resonant  photopumping  has  now  been  demonstrated  on  3  different  pulsed-power 
machines,  ranging  in  peak  current  from  0.4  MA  (LION,  Cornell),  to  1  MA  (Gamble  II, 
NRL),  to  10  MA  (Saturn,  Sandia  National  Laboratories).  The  Sandia  experiment  also 
produced  population  inversion. 

Ultra-intense,  femtosecond-pulsed  laboratory  lasers  ranging  from  the  ultraviolet  to 
the  infrared  represent  an  exciting  new  tool  for  investigations  varying  from  basic  atomic 
and  plasma  physics  to  specific  mission-oriented  goals.  We  have  continued  our  assessment 
of  the  potentialities  of  these  devices  with  respect  to  their  production  of  gain  at  ultrashort 
wavelengths  and  also  for  excitation  of  superintense  incoherent  x-ray  sources.  Section  2 
considers  the  interaction  of  a  subpicosecond  laser  pulse  with  a  planar  aluminum  target, 
using  a  model  which  includes  solution  of  the  Helmholtz  wave  equation  as  well  as 
self-consistent  detailed  configuration  atomic  physics  coupled  to  radiation  hydrodynamics. 
Among  other  results,  the  temperature  and  density  profiles  midway  from  the  leading  edge 
of  the  blowoff  plasma  are  found  to  support  population  inversions  for  a  number  of  upper 
states.  A  remaining  question  is  whether  these  inversions  persist  for  timescales  of  interest 
for  extracting  gain.  In  many  laser-produced  plasmas,  the  most  intense  radiation  emanates 


from  the  strongest  resonance  lines  which  by  nature  of  their  large  oscillator  strengths  are 
most  likely  to  be  optically  thick.  Escape  of  this  radiation  is  a  diffusion-like  process 
which  can  affect  the  x-ray  pulsewidth  and  thus  the  intensity.  In  section  3,  the  basic 
time-dependent  transfer  model  which  was  described  in  last  year’s  report  is  applied  to  the 
He-like  a  line  which  is  often  the  most  intense  K  shell  radiator.  Wide-ranging  results  for 
the  pulsewidth  as  well  as  analytic  models  including  atomic  number  scaling  enable  accurate 
assessments  to  be  made  of  this  effect  in  virtually  any  accessible  laboratory  experiment. 

Various  laser  absorption  mechanisms  can  result  in  the  production  of  non-Maxwellian 
electron  distributions  which  in  turn  affect  the  atomic  state  populations  by  altering  the 
relevant  rates.  This  problem  is  considered  in  detail  in  Section  4  for  neonlike  lasers  (the 
most  robust  of  the  short  wavelength  systems).  It  is  found  that  generation  of  ion-acoustic 
turbulence  at  high  laser  intensities  can  in  principle  increase  the  gain  of  3s-3p  transitions 
and  improve  the  performance  of  these  lasers. 


Enhancement  of  Ionization  in  Neon  X-ray  Laser  Targets  Radiatively 
Driven  by  a  Sodium-Potassium  Z  Pinch 

T.  Nash,  R.  Spielman,  and  T.  Tanaka 
Sandia  National  Laboratories 

J.  Apruzese  and  R.  W.  Clark 
Naval  Research  Laboratory 


Abstract 

We  present  spatially  and  temporally  resolved  XUV  spectra  from  a 
neon  gas  cell  lasant  target  that  is  photopumped  by  a  high  power 
sodium  z  pinch  2  cm  away  radiating  10  TW.  In  this  sodium/neon 
photoresonant  laser  scheme  the  2  to  1  transition  of  He-like  Na  in  the 
pinch  nominally  radiates  200  GW  and  is  resonant  with  the  4  to  1 
transition  in  He-like  neon.  We  have  measured  the  ionization  balance 
of  the  neon  gas  cell  as  a  function  of  both  space  and  time  and  found 
that  with  a  pure  sodium  pinch  the  neon  L  shell  is  not  completely 
ionized.  By  introducing  25%  potassium  into  the  z  pinch  load  we  have 
increased  the  amount  of  He-like  neon  in  the  target  without 
appreciably  affecting  the  sodium  pump  power.  The  spatial  profile  of 
neon  ionization  states  is  consistent  with  a  radiatively  driven  shock 
propagating  through  and  heating  the  neon  gas  cell  target. 


Manuscript  approved  March  12,  1993. 
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I.  Introduction 


X-ray  lasers  demonstrated  to  date  have  been  pumped  by  collisional 
excitation  and  recombination.  1-10  Although  resonant 
photopumping  could  potentially  produce  an  efficient  x-ray  laser  11- 
20,the  shortest  wavelength  at  which  photopumped  gain  has  been 
measured  is  2163  A  in  Be-like  C.  21 

Using  elliptical  crystal  spectroscopy,  resonant  photopumping  of  the 
n=4  level  of  He-like  neon  by  a  sodium  z  pinch  on  the  pulsed  power 
driver  Saturn  has  been  demonstrated.  22  The  200  GW  pump  is 
sufficient  to  produce  a  measureable  4  to  3  inversion  in  He-like  neon. 
Details  are  presented  in  Reference  22.  The  next  logical  step  in  this  x- 
ray  laser  research  is  to  measure  the  strongest  predicted  laser 
transition,  the  4f-3d  singlet  line  at  231  A. 

The  atomic  energy  level  diagram  including  the  coincident  lines  for 
the  sodium/neon  x-ray  laser  scheme  is  shown  in  figure  1.  The 
resonant  photopump  of  typically  200  GW  in  the  He-like  sodium  2  to 
1  transition  is  provided  by  the  z  pinch  load  of  the  Saturn  accelerator. 
Saturn  also  emits  up  to  10  TW  total  radiation  which  photoionizes  the 
neon  L  shell  in  the  neon  lasant  target,  and  is  predicted  to  drive  a 
shock  wave  from  the  side  window  of  the  gas  cell  into  the  neon.  23 
The  4p  singlet  level  of  He-like  neon  is  resonantly  photopumped. 
Electron  and  ion  collisions  distribute  the  population  and  hence  the  4 
to  3  inversion  over  the  n=4  sublevels  and  the  largest  gain  is 
predicted  for  the  4f  to  3d  singlet  transition  at  231  A.  13,19 

The  experimental  arrangement  is  shown  in  figure  2.  Saturn  delivers 
approximately  10  MA  with  a  60  nsec  rise  time  to  a  sodium  z  pinch 
load.  The  pinch  radiates  up  to  10  TW  with  200  GW  in  the  resonant 
photopump.  The  pulsewidth  of  the  resonant  photopump  is  20  to  30 
nsec.  A  neon  gas  cell  is  placed  with  its  side  window  18  mm  from  the 
pinch  axis,  just  outside  of  the  pinch  current  return  posts.  This 
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He-like  Na 
(pump) 


He-like  Ne 
(lasant) 


1.  Atomic  energy  level  diagram  for  Sodium/Neon  photoresonant  x- 
ray  laser  scheme. 


2.  Experimental  Layout  of  Z-pinch,  target,  and  diagnostic. 
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assures  that  the  neon  is  heated  by  radiation  rather  than  current.  A 
neon  fill  pressure  of  10  torr  is  confined  by  a  4000  A  Lexan  side 
window  and  an  end  window  that  is  either  2000  A  Lexan  or  5000  A 
polystyrene.  The  pumping  radiation  is  incident  on  and  ionizes  the 
side  window,  and  also  heats  the  neon  gas  target.  XUV  radiation  is 
measured  along  an  axial  line  of  sight  passing  through  the  end 
window  with  a  spatially  resolved  10  nsec  gated  MCPIGS 
(microchannel  plate  (MCP)  intensified  one  meter  grazing  incidence 

spectrometer). 2  4 


II.  Instrumentation 

We  measure  the  XUV  spectra  from  the  neon  gas  cell  with  a  spatially 
resolved  single  time  gated  MCPIGS.  The  MCP  is  curved  to  the  one 
meter  diameter  of  the  Rowland  circle.  The  instrumental  acceptance 
is  f/40  up  to  a  3  meter  radius  2  degree  grazing  incidence  spherical 
gathering  mirror. 

A  6  mm  wide  by  150  mm  long  stripline  on  the  MCP  is  gated  by  a  10 
nsec  900  volt  pulse  to  provide  a  single  10  nsec  resolved  time  frame. 
Typically  the  gate  is  centered  on  the  peak  of  the  z  pinch  x-ray 
emission.  On  the  short  6  mm  width  of  the  strip  we  spatially  resolve 
with  an  imaging  cylindrical  mirror  at  5  degrees  grazing  incidence. 
With  spatial  magnification  of  1/2.7  the  spatial  resolution  at  the 
target  is  better  than  1  mm.  However  due  to  the  large  150  mm  length 
of  the  spectral  focus,  the  spatial  focus  degrades  rapidly  away  from 
the  spatial  focus  wavelength.  At  spectral  locations  15  A  from  the 
spatial  focus,  the  spatial  resolution  degrades  to  worse  than  1  mm. 

This  limitation  in  the  spatial  imaging  of  the  MCPIGS  will  be 
addressed  in  the  future  by  using  a  flat  field  variable  line  spaced 
grating  spectrometer  that  maintains  spatial  focus  along  the  entire 
spectral  range.  On  the  MCPIGS  instrument  we  use  two  vacuum 
spacers,  one  for  spatial  focus  at  230  A  and  one  for  spatial  focus  at  80 
A,  assuming  the  use  of  a  typical  1200  1/mm  grating. 
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For  the  spatial  focus  at  230  A  we  use  a  1200  1/mm  grating  blazed  at 
5  degrees,  and  a  1000  A  AI  filter  to  remove  higher  orders.  For 
spatial  focus  at  80  A  we  use  no  filter  and  a  1200  1/mm  grating 
blazed  at  2  degrees.  The  entrance  slit  width  is  typically  20  microns 
and  the  grating  is  half  masked. 

A  spherical  spectral  gathering  mirror  sets  the  field  of  view  as  1  mm 
wide  at  the  gas  cell,  the  other  dimension  of  the  gas  cell  being  imaged 
by  the  cylindrical  mirror.  A  close  to  close  fast  valve  with  3  mm  by  3 
cm  aperture  open  for  500  microseconds  protects  the  mirrors  from 
shot  debris. 

The  microchannel  plate  is  backed  by  a  fiber  optic  faceplate  coated 
with  Pll  phosphor.  The  spectral  data  is  recorded  on  Kodak  type 
2484  film,  then  digitized  and  computer  processed. 


III.  Results 

Our  initial  and  unsuccessful  attempt  at  measuring  the  231  A  line  is 
depicted  in  figure  3.  When  viewed  through  a  2000  A  thick  lexan  end 
window  the  spectrum  is  plagued  with  oxygen  absorption  lines,  one  of 
which  falls  precisely  on  the  location  of  the  laser  line.  We  then 
eliminated  these  oxygen  lines  by  using  a  5000  A  oxygen-free 
polystyrene  end  window.  With  this  end  window  the  oxygen 
absorption  lines  disappear  and  the  background  level,  likely  due  to 
oxygen  lines  as  well  as  continuum  emission,  is  reduced.  Estimates  of 
instrumental  sensitivity  and  fluorescent  yield  in  the  laser  line 
indicate  that  we  are  unlikely  to  detect  the  231  A  transition  unless  it 
is  lasing  with  a  gain  greater  than  1  cm'  1 . 

Measurements  of  ionization  balance  are  made  by  recording  spectra  in 
the  70  to  120  A  window.  This  region  includes  bright  3  to  2 
transitions  in  Be-like,  Li-like,  and  He-like  neon.  With  a  pure  sodium 
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polystyrene  end  window  eliminates  oxygen  absorption. 


z  pinch  pump  the  lines  of  Be-like  and  Li-like  neon  are  about  equal  as 
depicted  in  the  bottom  trace  of  figure  4.  One  reason  we  are  not 
detecting  the  Lser  line  could  be  a  paucity  of  He-like  neon.  Complete 
photoionization  of  the  neon  L  shell  requires  copious  radiation  at  240 
to  450  eV  energies  ,  and  a  sodium  plasma  is  lacking  in  transitions 
near  these  energies.  In  order  to  increase  the  photoionization  pump 
we  introduce  potassium  into  the  z  pinch  with  the  idea  that  potassium 
L  radiation  could  help  photoionize  the  neon  L  shell  and  produce  more 
He-like  neon.  With  25%  potassium  in  the  driving  pinch  the  ratio  of 
Li-like  to  Be-like  neon  line  intensities  increases  a  factor  of  2 
indicating  a  higher  ionization  level  and  thus  more  He-like  neon  is 
present  when  potassium  is  used  in  the  pinch. 

The  spectra  of  both  pure  sodium  and  sodium-potassium  pinches 
were  measured  with  a  grazing  incidence  spectrometer.  This  data  is 
shown  in  figure  5.  The  spectral  distributions  are  based  on  an 
estimate  of  instrumental  response.  The  signal  at  all  subkilovolt 
energies  increases  a  factor  of  3  when  25%  of  the  pinch  mass  is 
potassium.  This  increases  the  ability  of  the  pumping  z  pinch  to 
photoionize  the  neon  L  shell,  as  corroborated  by  the  data  of  figure  4. 

The  effect  of  different  amounts  of  potassium  introduced  into  the 
sodium  pinch  was  scanned  on  a  shot  to  shot  basis.  The  results  are 
shown  in  figure  6.  The  total  number  of  wires  in  the  z  pinch  load 
array  is  16.  The  power  in  the  250  to  450  eV  window  is  measured  by 
a  titanium-filtered  carbon  x-ray  diode.  This  diode  shows  the  power 
in  this  window  increasing  a  factor  of  3  in  going  from  0  to  25% 
potassium  in  the  load,  and  is  in  agreement  with  the  spectroscopic 
data  of  figure  5.  The  1 1  A  resonant  photopump  line  power  and  total 
sodium  K  shell  power,  measured  by  x-ray  diodes,  bolometers,  and 
photoconductive  detectors,  do  not  signifigantly  decrease  up  to  a  25% 
fraction  (4  wires)  of  potassium.  However  at  50%  potassium  (  8 
wires)  the  sodium  K  shell  and  resonant  photopump  power  fall  off 
dramatically.  The  data  of  figure  6  shows  that  25%  potassium  in  the 
sodium  load  optimizes  the  ability  of  the  z  pinch  to  both  photoionize 
the  neon  L  shell  and  resonantly  photopump  the  He-like  neon. 
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A  small  amount  of  potassium  impurity  greatly  increases  the  soft 
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6.  Four  potassium  wires  maximizes  the  250-450  eV  power  without 
appreciably  affecting  the  pump  power. 


IV.  Spatial  Ionization  Profiles  and  Interpretation 


The  spatial  distributions  of  lines  of  He-like,  Li-like,  and  Be-like  neon, 
and  H-like  carbon,  recorded  with  a  10  nsec  gate  centered  on  the 
radiation  power  pulse  are  shown  in  figure  7.  The  lines  are  the  78.3 
A  3d-2p  triplet  transition  for  He-like  neon,  the  98.2  A  3d-2p 
transition  for  Li-like  neon,  the  cluster  of  3d-2p  triplet  transitions 
around  110  A  for  Be-like  neon,  and  Balmer  alpha  at  182  A  for  H-like 
carbon.  Be-like  and  Li-like  neon  are  seen  to  bum  out  at  the  center  of 
the  gas  cell  where  He-like  emissions  peak.  The  neon  is  likely  heated 
by  a  radiatively  driven  shock  wave  that  moves  from  the  side 
window  across  the  gas  cell.  The  H-like  carbon  in  figure  7  comes  from 
the  side  window  and  its  location  of  intensity  fall-off  indicates  the 
position  of  the  shock.  The  shock  is  estimated  to  travel  1  mm  during 
the  10  nsec  data  acquisition  gate.  The  shock  evidently  heats  a 
localized  region  of  the  neon,  burning  through  the  L  shell  and 
enhancing  He-like  neon  in  concert  with  direct  photoionization. 

As  mentioned  above,  the  He-like  neon  line  measured  is  the  78.3  A  3d 
to  2p  triplet  transition  array.  This  is  the  only  strong  XUV  line  from 
He-like  neon  we  have  detected.  Modelling  indicates  that  this  line 
should  be  stronger  than  other  He-like  neon  XUV  lines  because  the 
lower  2p  triplet  level  of  this  transition  has  a  slow  decay  time  and  a 
large  population.  We  have  also  detected  the  3p-2s  triplet  transition 
at  74.4  A  in  He-like  neon  at  a  level  barely  above  background  noise. 

The  small  measured  intensity  of  the  3d-2p  He-like  line  points  out  the 
difficulty  of  measuring  the  candidate  laser  transition  which  is 
predicted  to  be  about  10  times  weaker  in  fluorescence. 

The  optimum  fill  pressure  for  detecting  He-like  neon  lines  in  the  XUV 
has  been  found  to  be  10  torr.  In  figure  8  the  upper  trace  shows  the 
78.3  A  He-like  neon  line  at  almost  1/2  the  intensity  of  the  Li-like  4d 
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7.  At  the  shock  front  the  He-like  line  peaks  where  the  Li-like  and  Be¬ 
like  lines  burn  out. 
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8.  The  He-like  line  is  strong  at  the  shock  front,  but  increasing 
pressure  causes  recombination. 


to  2p  line,  measured  at  the  shock  front  with  a  10  torr  fill  pressure. 
At  30  torr  fill  pressure  the  78.3  A  line  intensity  increases  only 
slightly  and  is  overwhelmed  by  surrounding  lines  from  Li-like  and 
Be-like  neon.  Modeling  indicates  that  the  higher  fill  pressure  causes 
recombinations  back  into  lower  ionization  states,  whereas  the 
ionization  continues  to  be  dominated  by  the  fixed  photon  source. 

The  evolution  of  the  neon  gas  cell  target  has  been  calculated  for  a 
variety  of  assumed  pump  pinch  spectra  using  the  one-dimensional 
planar  radiation-hydrodynamics  model  described  in  Reference  25. 

In  these  calculations,  the  assumed  initial  neon  fill  pressure  has  also 
been  varied.  The  neon  atomic  model  includes  52  spectral  lines  in  Ne 
VIII  through  X  and  29  continuum  bound-free  edges,  all  self- 
consistently  coupled  to  the  energy  equatiion  via  multifrequency 
radiation  trasnport.  Some  of  the  initial  results  of  these  calculations 
have  been  described  in  Reference  23. 

Varying  both  the  assumed  pump  radiation  and  neon  fill  pressure 
obviously  leads  to  differences  in  the  calculated  evolution  of  the 
target  which  aid  in  interpreting  the  experimental  results.  However, 
the  basic  hydrodynamic  behavior  of  the  target  is  a  common  thread 
which  runs  through  all  of  the  numerical  solutions.  The  optical  depth 
of  the  cold  window  is  of  order  unity  to  the  soft  pinch  x-rays  in  the 
100  to  300  eV  spectral  region.  As  ionization  proceeds,  the  window 
plasma  becomes  increasingly  transparent  to  the  pinch  radiation, 
nonetheless  absorbing  enough  energy  to  strip  the  carbon  to  the  K 
shell,  and  behaving  like  a  laser-driven  exploding  foil.26  In  accord 
with  the  similarity  solutions  of  Reference  26,  the  present  numerical 
results  indicate  window  expansion  velocities  about  50%  greater  than 
the  local  isothermal  sound  speed  (3  x  10^  cm  sec"*  at  15  eV  to 
which  the  window  is  heated).  The  window  plasma  pushes  the  much 
lower  density  neon  1.5  to  2.0  mm  back  from  its  original  position 
during  the  approximate  30  nsec  x-ray  heating  to  the  peak  of  the 
pump  pulse.  This  creates  a  compressed,  heated  neon  region  about  1 
mm  thick  of  35  eV  temperature  contiguous  to  the  exploded  side 
window  plasma.  The  remainder  of  the  6  mm  neon  column  remains 
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quiescent,  heated  purely  by  x-rays  to  temperatures  of  10-15  eV. 

The  existence  of  a  30  eV  region  has  already  been  suggested  by 
experimental  measurements  of  the  slope  of  the  neon  recombination 
continuum. 22  Besides  demonstrating  the  enhancement  of  ionization 
due  to  the  potassium  radiation,  the  present  experimental  data  also 
verify  the  existence  of  the  predicted  shock-heated  neon  and  window 
explosion,  as  now  discussed. 

One  of  the  most  revealing  and  critical  areas  for  comparison  of  theory 
and  experiment  is  the  predicted  versus  measured  spatial  profile  of 
the  78.3  A  triplet  array  in  He-like  neon.  The  intensity  of  this 
transition  is  proportional  to  the  local  He-like  ground  state  density  as 
well  as  the  resonant  photopumping  rate  to  n=4.  No  collisional 
excitation  from  the  He-like  ground  state  occurrs  at  the  relatively  low 
temperatures  (<35  eV)  of  the  target,  whose  high  ionization  state  is 
principally  produced  by  photons  from  the  Z  pinch. 

In  Figure  9  the  calculated  n=3  triplet  population  density  at  the  peak 
of  the  x-ray  pulse  of  a  Na/K  pinch  is  plotted  versus  position  in  the 
gas  cell.  In  this  calculation  the  measured  spectral  shape  of  the  Na/K 
soft  x-ray  spectrum  was  used  to  drive  the  gas  cell  target.  The 
temporal  dependence  of  the  pulse  was  assumed  to  be  Gaussian  with 
a  FWHM  of  35  nsec.  Note  the  excellent  qualitative  correspondence 
between  Figure  9  and  the  He-like  trace  in  Figure  7.  Both  the 
measured  and  calculated  peaks  occur  3. 5 -4.0  mm  to  the  rear  of  the 
original  side  retaining  window.  The  calculation  indicates  that 
explosion  of  the  window  has  pushed  the  carbon-neon  interface  back 
1.8  mm  from  its  original  position.  The  spatially  resolved 
experimantal  data  of  Figure  7  also  indicates  shock  positions  of  1-2 
mm(some  blurring  occurs  from  integrating  over  the  10  nsec  gate). 

The  basic  hydrodynamic  processes  occuring  in  these  x-ray  driven  gas 
cell  targets  therefore  appear  to  be  well  understood  and 
characterized. 

Since  the  ground  state  ionization  thresholds  of  B-like,  Be-like,  and  Li- 
like  neon  lie  at  158,  207,  and  239  eV,  respectively,  it  is  clear  that  the 
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9.  The  calculated  He-like  n=3  triplet  population  is  plotted  versus 
position  in  the  gas  cell  at  the  peak  of  the  x-ray  pulse  from  the 
driving  sodium/potassium  pinch.  The  initial  fill  pressure  is  5  torr. 
The  original  position  of  the  gas  cell  side  window  is  at  0. 


large  increase  in  spectral  output  in  this  region  due  to  the 
introduction  of  potassium  will  enhance  the  ionization  of  the  neon. 
The  He-like  neon  fraction  is  calculated  to  vary  from  55%  for  a  pure 
Na  pump  with  10  torr  fill  pressure  to  80%  for  a  Na/K  pump  at  5  torr 
fill. 


V.  Conclusions 

We  have  measured  the  ionization  balance  across  a  neon  gas  cell 
irradiated  by  a  10  TW  z  pinch.  The  measurements  are  consistent 
with  a  radiatively  driven  shock  wave  propagating  through  the  neon. 
We  have  increased  the  overall  level  of  ionization  and  the  amount  of 
He-like  neon  in  the  gas  cell  by  introducing  25%  potassium  into  the  z 
pinch  load.  This  increases  the  ionizing  soft  x-ray  radiation  a  factor  of 
3.  Oxygen  absorption  lines  have  also  been  eliminated  from  the 
spectrum  by  using  an  oxygen-free  5000  A  polystyrene  end  window. 
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Ultra  Short  Pulse  Laser  Produced  Aluminum  Plasma 
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ABSTRACT 

The  interaction  of  a  subpicosecond  laser  pulse  with  a  planar  aluminum  target  is  investigated. 
The  interaction  is  described  using  a  1-D  time-dependent  non-LTE  radiat.jn  hydrodynamics 
model  in  conjunction  with  a  Helmholtz  wave  equation  characterizing  the  laser-target  coupling.  In 
addition,  a  fraction  of  the  absorbed  laser  energy  is  partitioned  into  energetic  electrons  that  are 

transported  and  deposited  into  the  cold  un-ionized  portion  of  the  target.  The  radiative 
characteristics  of  the  Ka  and  the  K-shell  line  spectra  are  presented  and  discussed.  The 

feasibility  of  creating  conditions  conducive  to  population  inversion  in  selected  transitions  in 
lithiumlike  aluminum  is  investigated. 


INTRODUCTION 

The  interaction  of  a  subpicosecond  pulsed  laser  with  a  solid  target  produces  a  hot  solid 
density  plasma  with  a  radiation  signature  that  is  characterized  by  an  ultra  short  burst  of 
energetic  x-rays.  The  duration  of  the  emitted  radii,  jon  pulse  is  on  the  order  of  picoseconds  and 
represents  a  potentially  useful  source  of  x-rays  of  interest  for  a  variety  of  applications  involving: 
ultra  fast  optics,  picosecond  spectroscooy,  inner  shell  photo-pumping  and  recombination  x-ray 
lasers,  and  chirped  pulsed  x-rays.  However,  before  this  new  tool  can  be  utilized  we  must  first 
understand  and  characterize  the  properties  and  behavior  of  the  plasma  that  generates  the  x- 
rays.  This  raises  a  number  of  issues,  the  first  and  most  important  involves  describing  how  the 
laser  energy  couples  to  the  cold  solid  target.  Conventional  wisdom  suggests  that  for  laser 
intensities  less  than  about  101®  watts/cm2  the  dominant  absorption  mechanisms  are  thought  to 
be  multiphoton  ionization  initially  followed  and  dominated  by  a  combination  of  inverse 
bremsstrahlung  and  resonance  absorption.  Above  10^®  watts/cm2  absorption  is  more  likely 
dominated  by  a  variety  of  parametric  plasma  processes.  Unfortunately,  in  general,  there  is  no 
unanimity  on  what  is  the  dominant  absorption  process  and  so  it  remains  an  open  issue.  A 
series  of  well  thought  out  and  well  diagnosed  experiments  would  provide  a  useful  data  base 
and  starting  point  for  addressing  the  absorption  issue. 
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In  this  paper  we  treat  the  laser  plasma  coupling  for  intensities  less  than  1018  watts/cm2 
and  for  S  polarization  by  means  of  a  Helmholtz  wave  equation  with  an  enhanced  collision 
frequency.  An  enhanced  collision  frequency  is  employed  in  order  to  guarantee  sufficient 
absorption  and  subsequent  plasma  heating  to  reproduce  the  observed  spectra.  In  addition,  in 
order  to  explain  some  of  the  observed  features  associated  with  the  inner  shell  x-ray  spectra,  we 
have  tacitly  assumed  that  a  fraction  of  the  absorbed  energy  is  partitioned  into  energetic 
electrons.  These  are  transported  and  deposited  in  the  un-ionized  cold  aluminum,  producing 
local  shock  heating  and  inner  shell  vacancies.  The  K-shell  vacancies  are  filled  rapidly  by  either 
Auger  decay  or  resonance  fluorescence. 

Numerical  simulations  are  carried  out  for  several  KrF  laser  intensities  ranging  from  1016 
to  3x1 01 7  watts/cm2  with  pulse  lengths  of  650  femtoseconds  normally  incident  on  a  2 p  thick 
planar  aluminum  slab.  Comparisons  are  made  for  two  laser  intensities,  two  collision 
frequencies,  and  two  beam  energies. 

MODEL 

The  calculations  are  based  on  a  time-dependent  1-D  collisional  radiative  radiation 
hydrodynamics  model  self-consistently  coupled  to  a  Helmholtz  wave  equation,  describing  the 
absorption  of  the  laser  pulse  which  is  assumed  normally  incident  on  the  target.  The  collisional 
radiative  model  contains  all  the  ground  states  and  selected  excited  levels  distributed  throughout 
the  manifold  of  ionization  stages.  In  particular,  the  K-shell  manifold  contains  excited  states  up  to 
and  including  the  n-10  level  and  contains  both  singlet  and  triplet  states  in  the  heliumlike  stage. 
The  lithiumlike  ionization  stage  also  includes  all  levels  up  to  n-10.  The  remaining  ionization 
stages  contain  ground  states  and  a  variety  of  exerted  states  important  to  the  radiative 
energetics  and  diagnostics  of  the  plasma.  The  ionization  dynamics  model  is  self-consistently 
coupled  to  a  multi-cell  probability  of  escape  radiation  transport  scheme.  The  spectral  line  profile 
functions  are  represented  by  Voigt  functions.  A  detailed  discussion  of  the  ionization  dynamics 
and  atomic  processes,  including  the  inner  shell  processes,  are  given  in  references  1  -3.  The 
absorption  of  the  normally  incident  laser  pulse  on  the  target  is  calculated  by  solving  the 
Helmholtz  wave  equation  for  a  multi-layered  medium.  In  this  way  we  self-consistently  allow  for 
the  interaction  of  the  incident,  reflected  and  transmitted  waves  in  each  numerical  zone  in  the 

medium.  The  resistivity  was  calculated  using  the  collision  frequency  model  of  Cauble  and 
Rozmus4,  which  is  represented  by  vo-  A  fraction  of  the  absorbed  energy,  in  our  case  30%,  was 

partitioned  into  energetic  electrons.  These  electrons  were  sic  to  down  and  deposited  their 
energy  according  to  the  Bethe-Bloch  approximation.  The  situations  were  performed  for 
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energetic  electrons  consisting  of  either  a  10  KeV  or  4  KeV  beam  energy  and  for  collision 
frequencies  10X  and  100X  ^respectively. 


RESULTS 

The  rate  at  which  a  Gaussian  profile  laser  pulse  of  FWHM  650  fsec  248  nm  KrF  laser 
pulse  of  10^6  watts/cm2  is  deposited  into  a  2 p  thick  aluminum  target  is  shown  in  Fig.  1  as  a 
function  of  position  at  0.4  psec.  The  laser  pulse  is  assumed  incident  on  the  target  from  right  to 
left.  At  this  time  about  14%  of  the  laser  power  is  deposited  in  the  target.  Also  shown  in  Fig.  1 
are  the  mass  density  and  the  ratio  of  the  local  electric  field,  averaged  over  an  optical  cycle, 
<E>,  to  the  amplitude  of  the  incident  electric  field,  Eg,  as  functions  of  position.  The  location  of 
the  critical  surface,  i.e.  where  co-cop  is  represented  as  Rcriticai*  There  is  an  enhancement  of 
<EWEo  at  about  2.08[i  that  represents  the  constructive  interference  between  the  incident  and 
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Fig.  1.  Laser  energy  deposition  at  0.40  picoseconds. 

The  logarithm  of  the  rate  of  energy  deposition  into  the 
target  BEi^^r/di  is  plotted,  and  the  location  of  the  critical 
surface  R is  shown.  The  local  electric  field,  averaged 
over  an  optical  cycle,  <  E  >  is  shown  normalised  to  the  am¬ 
plitude  of  the  incident  field  E ..  The  mass  density  p  is  shown 
in  a  stepwise  manner  to  indicate  the  spatial  resolution  of  the 
simulation  in  the  deposition  region. 

Displacement  (p) 


Fig.  2.  Hydrodynamic  profile*  at  0.40  picosecond*. 

Spatial  profiles  of  mass  density  and  temperature  are  plot¬ 
ted  near  the  peat  of  the  User  pulse.  Two  effective  collision 
frequencies  are  compared  for  a  laser  intensity  of  10“  IV/cm’: 
10f.  and  100o.. 


reflected  waves.  The  mass  density,  log(p),  is  plotted  in  "  staircase"  mode  to  emphasize  the 
spatial  resolution  of  the  numerical  grid.  At  this  time  only  a  very  weak  shock  wave  due  to  the 
deposition  appears  at  about  1 .86p.  Figure  2  shows  the  hydrodynamic  profiles  at  0.4  psec  for  1 0 
vg  and  IOOvq  as  a  function  of  distance.  These  profiles  include  the  effects  of  the  deposition  due 
to  an  energetic  beam  of  10  KeV  electrons  containing  30%  of  the  absorbed  energy.  Since  there 
is  little  expansion  of  the  leading  edge  of  the  blowoff  plasma  at  this  early  time  the  density  profiles 
remain  essentially  the  same.  On  the  other  hand,  the  temperature  profiles  are  significantly 
different.  As  a  result  of  the  enhanced  IOOvq  collisional  and  beam  heating  a  higher  temperature 
profile  is  produced  both  in  the  laser  deposition  region  and  deeper  into  the  target  than  in  the  1 0x 
vg  case.  Also,  as  a  consequence  of  the  beam  deposition  there  is  an  additional  local  heating 
causing  a  second  weak  shock  wave  deep  in  the  cold  aluminum.  The  temperature  differential 
between  these  two  collisional  cases  is  simply  the  result  of  distributing  more  of  the  absorbed 
energy  into  thermal  energy  in  the  IOOvq  case  than  the  10vq  case.  In  the  absence  of  an 
energetic  electron  beam,  the  temperature  profile  initially  would  be  dominated  mostly  by  thermal 
conduction  from  the  laser  deposition  region  and  exhibit  a  rapid  decay  into  the  cold  aluminum. 
These  differences  are  easily  observed  and  identified  in  the  x-ray  spectra. 

In  Figures  3  and  4  are  shown  the  spatially  integrated  and  time  resolved  intensity  of  the  x- 
ray  emission  spectra  for  the  lines  and  K-shell  lines  as  a  function  of  photon  energy  at  0.4 

psec  for  an  incident  laser  intensity  of  1016  watts/cm2  with  a  10  KeV  beam  for  10  and  IOOvq 

respectively.  The  lines  are  identified  and  labeled  on  the  Figures.  The  most  notable  and 
distinguishing  feature  on  these  Figures  is  the  series  of  Kq,  lines  generated  by  the  energetic 

electrons.  Another  obvious  feature  on  Fig.  4  is  the  enhancement  (  by  almost  an  order  of 
magnitude  )  of  the  free-bound  continuum  which  is  observable  over  the  entire  photon  range 
shown.  Only  a  brief  glimpse  of  the  free-bound  continuum  is  visible  around  1 .4  KeV  on  Fig.  3. 
The  differences  in  the  intensity  of  the  lines  are  due  to  the  enhanced  collisionality.  The 

primary  reason  for  this  is  the  change  in  the  plasmas'  reflectivity.  Although  the  beam  energy  is 

the  same  in  both  cases  the  collisionality  differs  by  a  factor  of  ten  between  the  two  cases  leading 
to  a  factor  of  3  increase  in  the  absorption  properties  of  the  plasma.  Spectroscopically,  the  Ka 

signatures  and  their  intensity  distribution  represent  a  useful  diagnostic  for  determining  the 
effects  of  the  energetic  electrons.  They  also  represent  a  substantial  fraction  of  the  radiated  x- 
ray  line  energy.  The  other  interesting  feature  appearing  on  these  Figures  is  the  presence  of 
some  of  the  K-shell  lines,  in  particular  the  hydrogen-like  lines.  For  the  IOOvq  case  the  spectra 
contains  lines  from  both  the  helium-  and  hydrogen-like  ionization  stages  with  the  heliumlike 
resonance  line  dominating  the  K-shell  spectra.  The  hydrogenlike  resonance  line  is  about  equal 
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(total)  spectrum  0.4000  psec 


(total)  spectrun  0.4000  psec 
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Fig.  S.  Emission  spectrum  at  0.40  pace  for  v  «=  lOi/.. 

Instantoneoos  emission  spectrum  near  the  p— V  of  the 
laser  pulse  showing  magnitudes  of  the  K,  and  K-sheO 
The  collision  frequency  has  been  enhanced  by  a  factor  of  10. 
Sdium-o  emission  and  Km  radiation  from  the  lower 
states  of  aluminum  dominate  the  spectrum. 


Energy  (keT) 


Fig.  4.  Emission  spectrum  at  0.40  psec  for  v  <=  I00*v 

Instantaneous  emission  spectrum  near  the  peak  of  the 
laser  pulse  with  the  collision  frequency  enhanced  by  n  fac¬ 
tor  of  100.  The  redaction  in  the  reflectivity  produces  s  hotter 
plasms,  evidenced  by  the  presence  of  H-fibe  mwm*  and  en¬ 
hanced  Km  radiation  from  aluminum  V-DL 


in  magnitude  to  the  Heliumlike  p  line  indicating  that  the  temperature  of  the  emitting  region  is 
insufficient  to  populate  the  hydrogenlike  levels  and  is  more  readily  conducive  to  populating  the 
upper  levels  of  the  heliumlike  states.  Both  are  significantly  less  than  the  heliumlike  resonance 
line  by  about  a  factor  of  50.  For  the  10vo  case  the  plasma  is  radiating  predominantly  from  the 
heliumlike  resonance  line.  Ultra  short  pulse  laser  aluminum  target  experiments  performed  at 
lO1^  watts/cm2  produce  an  emission  spectra  much  like  the  IOvq  case  where  hydrogenlike 
spectra  are  absent  and  only  the  heliumlike  and  Ka  lines  appear.  5 

Figures  5  and  6  represent  the  total  integrated  spatial  and  spectral  intensity  up  to  620 
femtoseconds  for  an  incident  laser  intensity  of  1017  watts/cm2  and  an  enhanced  collisionality  of 
IOOvq  for  electron  beam  energies  of  4  and  10  KeV,  respectively.  The  significant  spectral 
differences  between  these  two  cases  is  the  intensity  and  distribution  of  the  Ka  lines  and  the 

magnitude  of  the  underlying  continuum.  These  observations  could  be  used  as  a  diagnostic  to 
infer  plasma  conditions.  The  heliumlike  resonance  line  radiates  like  a  beacon  and  dominates 
the  K-shell  spectrum  in  both  cases.  However,  at  a  laser  intensity  of  10"17  watts/cm2  and  a  100 
fold  increase  in  the  collision  frequency  the  plasma  is  hot  enough  to  sustain  and  radiate 
significant  amounts  of  both  helium-  and  hydrogen-like  lines.  The  slight  difference  in  the 
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magnitude  of  the  underlying  continuum  between  the  two  cases  is  due  to  the  increased  free- 
bound  recombination  radiation  in  the  4  KeV  beam  case  due  to  the  cooler  temperature  in  the 

target  as  a  result  of  the  lower  beam  energy.  This  behavior  is  also  reflected  in  the  distribution 
and  intensity  of  the  Ka  lines. 


told  iategrtved  ipectrun  0.6290  pi«e 


Energy  (ke¥) 


Fig.  6.  Integrated  Eminian  spectrum  far  Eh..  =  4keV. 

Emission  spectrum  time-integrated  over  the  User  pulse 
for  4  keV  fast  electrons.  Km  radiation  from  the  lower  ionised 
states  of  aluminum  is  relatively  weak,  since  fast  electrons  are 
being  deposited  in  hot  plasma  near  the  critical  surface. 


total  integrated  sgectrisn  0.6200  pace 


Energy  (keV) 


Fig.  6.  Integrated  Emission  spectrum  for  Eh..  =  lOkeV . 

Emicakm  spectrum  time-integrated  oner  the  laser  pulse 
(or  10  keV  fast  electrons.  Km  radiation  from  the  lower  ionised 
states  of  aluminum  is  now  dominant,  since  many  energetic 
electrons  are  being  deposited  in  the  cool  interior  of  the  target. 


Finally,  we  have  investigated  the  possibility  of  plasma  conditions  being  generated  that 
are  conducive  to  creating  population  inversions  in  selected  levels  in  lithiumlike  aluminum.  For 
the  case  of  a  laser  intensity  of  3x10^  watts/cm^  with  a  100  fold  increase  in  collisionality  and  a 
beam  energy  of  10  KeV,  preliminary  results  indicate  that  the  temperature  and  density  profiles 
midway  from  the  leading  edge  of  the  blowoff  plasma  do  support  population  inversions  in  a 
number  of  upper  states.  As  an  example,  we  present  results  for  the  n=8  principal  quantum 
number  of  lithiumlike  aluminum  in  Frg.  7.  The  energy  levels  and  their  energies  relative  to  the 
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ground  state  are  shown  along  with  the  ratio  of  fractional  population  of  the  upper  to  lower  level, 
i.e.,  (fu/gu)/tyg|)  where  f u  |  and  gU||  represents  the  fractional  population  and  statistical  weight  of 
the  upper  and  lower  level,  respectively.  The  results  of  these  simulations  provide  encouraging 
data  supporting  the  idea  that  it  is  possible  to  create  inversions  over  timescales  of  interest. 
Whether  the  existence  of  population  inversions  can  successfully  lead  to  lasing  with  a  gain 
length  is  still  under  active  investigation. 
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Fig.  7.  Population  Inversions  in  Li-like  Aluminum. 

Ratio  of  Lithium-like  fractional  populations  f%gi/fig%  at 
0.48  picoseconds  for  x=1.82#i  (0.05p  behind  the  critical  sur¬ 
face).  Population  inversions  among  the  atomic  levels  (mainly 
TV  >  3)  are  produced  at  this  time.  Those  relative  to  the  N  =8 
level  are  shown  above. 
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ABSTRACT 


The  heliumlike  alpha  resonance  line  is  one  of  the  strongest  and  most  frequently  observed 
transitions  in  plasma  spectroscopy.  In  optically  thick  plasmas  of  moderate  density,  photons  emitted 
in  this  line  typically  undergo  many  absorptions  and  re-emissions  before  escaping.  In  the  present 
investigation  the  time  required  for  this  process  to  occur  in  aluminum  plasma  is  calculated  by 
numerical  solution  of  the  time-dependent  equation  of  radiative  transfer,  in  which  the  photon  transit 
time  between  interactions  is  taken  into  account.  An  analytic  model  is  also  developed  in  which 
the  collisional  quenching  probability  is  parameterized  in  terms  of  an  equivalent  two-level  atom. 
This  model  facilitates  interpretation  of  die  numerical  results,  and  permits  economical  estimates  of 
the  escape  time  for  different  plasma  elements  and  conditions.  One  practical  implication  of  these 
results  is  that  a  subpicosecond  x-ray  pulse  from  a  femtosecond  laser-produced  high  density  plasma 
can  be  broadened  by  multiple  scattering  if  it  propagates  into  a  moderate-density  preformed  plasma. 


L  INTRODUCTION 


The  ls2-ls2p1P1  principal  resonance  line  of  the  heliumlike  ion  (designated  He-a)  is  one  of 
the  most  frequently  observed  lines  in  plasma  spectroscopy.  It  has  been  detected  in  devices  ranging 
from  tokamaks  to  inertial  confinement  fusion  capsules  to  Z  pinches.  Most  recently,  this  line  has 
appeared  prominently  in  the  spectra  of  targets  irradiated  by  high  power  sub-picosecond  lasers. 

For  elements  with  atomic  number  Z  between  4  and  SO,  the  oscillator  strength  of  this  transition 
[1]  exceeds  0.5,  hence  it  is  usually  optically  thick  for  moderate  to  high  density  plasmas.  The 
associated  rapid  spontaneous  decay  rate  means  that  absorption  of  a  line  photon  is  often  followed 
by  re-emission  within  the  line  profile  (referred  to  as  scattering). Thus,  much  of  the  line  energy  that 
escapes  the  plasma  does  so  after  diffusion  in  both  position  and  frequency,  which  can  affect  the 
intensity  and  width  of  the  x-ray  pulse.  For  multinanosecond  sources  such  as  the  dense  Z  pinch,  this 
effect  is  expected  to  be  small  compared  to  the  inherent  lifetime  of  the  plasma.  However,  this  may 
not  be  the  case  for  the  very  short  lived  plasmas  produced  by  the  current  generation  of  powerful 
subpicosecond  lasers  [2-12].  Numerous  uses  of  such  plasmas  are  envisioned,  prominently  among 
them  the  production  of  ultrashort  x-ray  pulses  [2,4,5,12-16]. 

One  interesting  phenomenon  which  has  been  identified  and  quantified  to  some  degree  is  the 
effect  of  a  laser  prepulse  on  the  characteristics  of  the  x-rays  emitted  from  such  short-lived  plasmas. 
The  Los  Alamos  group  [7, 17]  and  others  [10]  report  that  the  presence  of  a  prepulse  has  a  strongly 
positive  effect  on  the  plasma  x-ray  yield.  This  occurs  whether  the  prepulse  is  of  nsec  [7,17]  or 
sub-psec  [10]  duration.  With  a  prepulse  present,  the  preformed  plasma  may  have  time  to  expand, 
and,  as  pointed  out  in  Refs.  14-16,  expansion  and  thermal  conduction  are  potentially  important 
cooling  mechanisms  in  these  plasmas.  The  greater  the  scale  length,  the  more  these  mechanisms 
tend  to  broaden  the  x-ray  pulse.  A  tradeoff  between  yield  and  pul  sc  width  thus  appears  to  be 
indicated. 

It  is  the  relatively  large  preformed  plasmas  producing  the  best  x-ray  yields  which  are  also  the 
most  likely  to  be  affected  by  pulse  broadening  inherent  in  the  multiple  scattering  of  strong  x-ray 
lines.  The  purpose  of  the  present  work  is  to  quantify  the  escape  time  of  the  He-a  line  radiation. 
Section  II  presents  the  numerical  and  analytic  models  of  radiation  transport  used  in  the  calculation. 
The  numerical  calculations  are  done  for  aluminum  which  has  been  a  frequent  experimental  target 
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for  sub-psec  lasers  [5,7,10,17],  and  from  which  the  He-a  line  has  been  strongly  detected.  In  Sec. 
m,  the  validity  of  the  analytic  model  is  assessed  by  direct  comparison  with  numerical  results. 
The  analytic  model  is  extendable  to  elements  other  than  aluminum.  Results  for  different  plasma 
conditions  are  also  presented  and  discussed. 

II.  RADIATION  TRANSPORT 

A.  Numerical  model 


The  time-dependent  transport  of  radiation  of  frequency  v  along  a  ray  in  the  positive  z  direction 
is  described  by  the  equation  of  transfer 


dlu  _  .  .  T  1  8IV 

dx  ~3v~kviv~  c~m' 


(1) 


In  eq.  (1),  J„  is  the  specific  intensity,  jv  and  are  the  emission  and  absorption  coefficients,  and  c 
is  the  velocity  of  light.  For  a  derivation  of  this  equation  see,  e.g..  Ref.  18.  Consider  the  spacetime 
point  (z0 ,  to  )•  Looking  in  the  negative  z  direction,  the  quantities  determining  J„(z0 ,  t0 )  are  jv  and 
kv  encountered  at  retarded  times  tr (z<j ,  to ,  z)  =  <o  —  .  The  effective  optical  depth  between 

zo  and  z  at  time  f0  is 


r*  o 

y(x0,t0,x)  =  J 


kv(tr)dx 


(2) 


and  the  solution  to  eq.  (1)  is 


M*o,<o)=/  jv[x,tT{x)]' exj>\-Tv{x0,tQ,x)\dx.  (3) 

In  eq.  (3),  z&  is  the  coordinate  of  the  plasma  boundary.  In  all  cases  to  be  considered,  it  is  assumed 
that  no  external  x-rays  are  incident  on  the  plasma.  Given  the  emission  and  absorption  coefficients 
from  previous  time  steps,  eq.  (3)  is  numerically  solved  on  a  grid  of  31  frequency  points  within 
the  line  profile  and  a  spatial  mesh  of  26-31  points  per  ray.  Linear  interpolation  in  both  space 
and  time  is  employed  for  jv  and  ku.  The  maximum  time  step  employed  is  10  fsec.  To  compute 
photoexcitation  and  escape  rates,  the  angle-averaged  intensity  and  net  radiative  flux  Fv  are 
needed.  This  requires  an  assumption  for  the  shape  and  dimensionality  of  the  medium.  The 
present  work  will  be  confined  to  one-dimensional  plane -parallel  geometry.  The  numerical  burden 
of  including  full  time  dependence  in  the  radiative  transfer  is  kept  at  an  acceptable  level  with  this 
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simple  choice  of  geometry.  Also,  the  optical  depth  presented  to  an  escaping  line  photon  can  be 
unambiguously  characterized  and  associated  with  a  calculated  escape  time  or  probability.  This 
will  facilitate  analysis  of  future  multidimensional  work  more  representative  of  actual  laboratory 
plasmas. 

The  distribution  of  the  rays  follows  Gaussian  quadrature  as  described  by  Chandrasekhar  in 
Ref.  19.  Each  ray  is  characterized  by  the  cosine  of  its  angle  (ji  =  cosd)  to  the  planar  normal.  In 
the  present  work  both  four  rays  and  two  rays  have  been  employed  for  comparison.  Their  specific 
/x’s  are  given  in  Table  HI  of  Ref.  19.  In  planar,  one-dimensional  geometry  the  quantities  J„  and 
F„  are  given  by 


J* 


(4) 


r+1 

Fv  =  2rr  /  Iv(fi)fidfi. 


(5) 


The  line  photon  escape  rate  is  the  frequency  integral  of  F„  at  the  plasma  boundary.The 
photoexcitation  rate  for  a  spectral  line  of  absorption  oscillator  strength  fiu  from  lower  level 
l  to  upper  level  u  is  given  by 


Wiuisec-1) 


(6) 


where  the  integral  is  taken  over  the  line  profile  <j>v  whose  normalization  is  J  <f>„dw  =  1.  Given  the 
photoexcitation  rates  as  well  as  the  appropriate  set  of  collisional  rates,  populations  of  the  various 
levels  and  ionic  stages  are  obtained  as  a  function  of  time  by  numerical  integration  of  the  atomic 
rate  equations  (see,  e.g.,  eq.  (1)  of  Ref.  20).  The  method  of  Young  and  Boris  [21]  is  employed  as 
it  is  especially  well  suited  to  the  time  integration  of  numerically  stiff  differential  equations. 


B.  Analytic  model 


The  two-level-atom  approximation  has  been  treated  extensively  in  the  astrophysical  literature 
(see,  e.g..  Ref.  22  and  references  therein).  It  is  therefore  a  natural  starting  point  for  an  analytic 
description  of  the  radiative  transfer  within  the  He-a  line.  Even  though  the  heliumlike  ion  is 
obviously  not  a  two-level  system,  the  same  processes  of  scattering  and  collisional  quenching  are 


32 


operative  in  determining  the  fate  of  the  emitted  line  photons.  Upon  absorption,  the  photon  will 
either  be  re-emitted  or  various  other  processes  will  transfer  the  ls2p*  Pi  electron  elsewhere.  For  the 
idealized  two-level-atom  case,  “elsewhere”  is  only  the  ground  state.  For  the  real  heliumlike  ion,  the 
electron  can  in  principle  end  up  almost  anywhere  within  the  same  ion  or  the  adjoining  hydrogenic 
or  lithiumlike  species.  One  cannot  obtain  a  probability  of  collisional  quenching  (designated  e)  by 
merely  adding  the  rates  to  all  other  possible  levels  from  ls2p:  Pi ,  since  in  some  cases  the  electron 
is  collisionally  shunted  right  back  to  the  original  upper  level  of  the  transition.  It  is  necessary  to 
subtract  (for  each  transition)  the  probability  that  the  electron  returns  to  the  ls2p*  Pi  level. 

If  the  electron  is  transferred  to  level  *  from  the  ls2p*  Pi  state  the  probability  of  return  to  the 
ls2pJPi  level  Pri  is  defined  as  the  rate  from  level  t  back  to  ls2p*Pi  divided  by  the  sum  of  all  of 
the  rates  out  of  level  i.  The  net  quenching  probability  e  is  then  given  by 


ZCjjl-Prt) 

Aul  +  ZCi(l-PriY 


(7) 


In  eq.  (7),  Ci  is  the  rate  from  ls2p*Pi  to  the  tth  level,  Aui  is  the  spontaneous  decay  rate  of  the 
resonance  transition,  and  the  sum  is  taken  over  all  levels  connected  to  ls2p*Pi . 

The  quenching  probability  e  has  been  computed  for  A1  using  the  atomic  data  described  in 
Ref.  20.  From  the  ls2s1S0  level  the  chance  of  return  to  the  ls2p1PJ  level  is  appreciable  (0.92  at 
an  electron  density  of  1021  cm-3  near  temperatures  of  200  eV).  At  these  same  conditions,  the  net 
probability  of  collisional  quenching  of  the  A1  He-a  line  per  absorption  (e)  is  0.029. 

Line  photons  may  also  be  eliminated  from  the  plasma  by  escape.  The  probability  Pe  that 
the  photon  escapes  a  medium  of  a  given  optical  depth  is  obtained  from  a  weighted  average  of 
the  exponential  escape  probability  at  each  frequency  within  the  line  profile.  The  weighting  factor 
is  the  strength  of  the  profile  itself  at  each  frequency.  Detailed  treatment  of  escape  probability 
calculations  may  be  found  in  Refs.  23-27.  For  stationary  media,  useful  analytic  forms  have  been 
obtained  for  Doppler  [23,24],  Voigt  [25],  and  Stark  [26,27]  line  profiles.  The  Doppler  shift  due  to 
bulk  motion  can  influence  Pe  in  some  instances.  For  media  where  this  is  important,  the  original 
work  of  Sobolev  [28]  has  recently  been  extended  by  Shestakov  and  Ed;r  [29]  to  include  cylindrical 
geometry. 

The  quantities  e  and  Pe  discussed  above,  along  with  the  spontaneous  decay  rate  of  the  upper 
level  Aui.  suffice  for  an  approximate  analytic  description  of  line  photon  escape  and  transport.  The 
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probability  that  a  photon  escapes  after  an  excitation  is  (1  -  e)Pe,  whereas  the  probability  that  it 
remains  within  the  plasma  is  5  =  (1  —  e)(l  —  ~4).  The  probability  that  the  photon  has  escaped  by 
the  N‘*  scattering  is  the  sum  of  the  probabilities  of  escape  after  initial  excitation  and  each  of  the 
subsequent  scatterings,  viz. 


Pn  =  (1  -  t)P* 


N 


i  +  £s" 

ns  1 


=  (1  -  e)Pe 


1-SN] 
1  -  5 


(8) 


The  ultimate  escape  probability  is  obtained  by  taking  the  limit  of  eq.  (8)  as  N— »  oo,  Pu  =  Pe(  ). 
The  time  required  for  N  scatterings,  i.e.,  the  time  associated  with  line  photon  probability  of  escape 
Pat,  including  the  initial  excitation  but  ignoring  the  photon  travel  time  between  absorptions,  is 
simply  the  level  lifetime  A~l  multiplied  by  N+l.  Thus,  Pn  given  by  eq.  (8)  is  equivalently  P(i), 
where  t  =  This  analytic  form  for  the  fraction  of  line  photons  which  have  escaped  as  a 

function  of  time  is  compared  in  the  next  section  with  numerical  results.  It  is  straightforward  to 
obtain  an  explicit  expression  for  the  time  for  fractional  escape  /e  (of  the  photons  which  ultimately 
escape).  The  condition  is  /*  =  1  —  SN,  which  leads  to 


and  in  the  limit  Pe,  e<l, 


<(/«)  =  Kl 
*</«)  =  Kl 


InS 


7n(l  -  /e) 

l-C^e  +  e) 


(9) 

(10) 


It  is  possible  to  construct  an  analytic  model  using  three  or  more  levels.  However,  as  pointed 
out  in  Ref.  30,  the  algebra  quickly  becomes  unmanageable.  Moreover,  as  seen  below,  the 
equivalent  two  level  approach  yields  respectable  accuracy  as  well  as  considerable  physical  insight. 
Note  also  that  a  single  Pe  is  employed  despite  the  fact  that  in  principle  it  varies  between  subsequent 
scatterings  when  the  photon  is  emitted  from  a  different  position  in  the  plasma.  However,  as  has 
been  shown  previously  [31,32],  the  line  photon  usually  remains  imprisoned  within  the  line  core 
and  stays  close  to  its  location  of  previous  emission,  escaping  in  one  long  flight  when  emitted 
in  the  line  wing.  Although  somewhat  oversimplified,  this  picture  explains  the  success  of  this 
approximation  in  both  the  present  and  previous  work.  As  evident  in  the  next  section,  some  spatial 
diffusion  of  the  excitation  does  occur. 


III.  RESULTS 
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A.  Equivalent  two-level -atom 


The  calculations  presented  in  this  section  have  two  principal  objectives.  First,  to  assess 
whether  under  the  chosen  conditions  transport  of  the  principal  resonance  line  of  the  heliumlike  ion 
can  be  understood  and  parameterized  as  if  the  ion  consisted  merely  of  a  ground  and  one  excited 
state.  Second,  to  establish  the  numerical  parameters  needed  to  insure  accuracy  in  the  computation. 

As  stated  above  in  Sec.  II B,  at  an  electron  density  of  1021  cm-8,  and  temperatures  around 
200  eV  the  net  probability  of  collisional  quenching  per  absorption  for  the  A1  He-a  line  is  0.029. 
A  two-level-atom  test  case  is  chosen  to  correspond  to  this  value.  At  this  density,  the  line  is  well 
described  by  a  Voigt  profile  whose  broadening  parameter  o=0.05.  A  planar  A1  plasma  of  width 
100  (im  would  have  about  90%  of  its  ions  in  the  heliumlike  stage  under  these  conditions.  The 
line  center  optical  depth  To  of  the  He-a  line,  measured  normally  from  the  midplane  to  plasma 
boundary,  is  90.  These  are  the  conditions  for  which  the  present  calculations  were  run;  obviously, 
no  extensive  atomic  data  file  is  required  for  this  simplified  case.  The  plasma  is  assumed  stationary 
(i.e.,  bulk  Doppler  motion  is  neglected  in  the  radiative  transfer).  To  examine  the  diffusion  and 
escape  of  line  photons  the  following  initial  excitation  conditions  are  assumed:  the  ions  are  entirely 
in  the  ground  state  except  for  the  spatially  innermost  12%  where  1%  of  the  ions  are  in  the  excited 
state,  see  Fig.  1.  The  clock  starts  at  t=0  and  the  line  photons  begin  escaping  and  migrating 
toward  the  plasma  boundary.  The  collisional  excitation  rate  is  set  to  zero  as  if  the  plasma  cools 
instantaneously.  In  this  fashion,  the  spatial  diffusion  and  escape  times  obtained  reflect  only  the 
effects  of  radiative  transfer. 

The  upper  level  fraction  as  a  function  of  space  and  time  is  illustrated  in  Figs.  1  and  2  for  5 
time  points  varying  from  0  to  7  psec,  for  both  4-ray  and  2-ray  Gaussian  quadrature.  Note  first  the 
excellent  agreement  between  the  n=2  and  n=4  cases,  strongly  suggesting  that  2  rays  are  adequate. 
Spatial  diffusion  of  the  upper  level  population  due  to  line  photon  absorption  and  re-emission  is 
clearly  evident.  Of  considerable  interest  is  the  rate  of  decay  of  the  upper  level  population  in  the 
central  zone.  If  the  plasma  were  optically  thin,  this  population  would  decrease  with  an  e-folding 
time  of  A'?  =36  fsec.  The  actual  lifetime  is  much  greater  in  the  optically  thick  medium  because  of 
re-excitation  of  the  level  by  line  photons,  as  is  evident  in  Figs.  1  and  2.  For  the  first  0.5  psec,  the 
effective  lifetime  of  the  upper  level  is  0.7  psec,  increasing  to  1.1  psec  for  the  final  3  psec  covered 
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FIO.  1.  A1  X3I  ls2p1PJ  population  fraction  is  shown  as  a  function  of  space  and  time  in  the 
two-level-atom  approximation  at  0,  0.5,  and  2.0  psec  following  initial  excitation.  The  plasma  is 
planar,  100  fim  thick,  with  an  assumed  temperature  of  200  eV  and  electron  density  10J1  cm-5. 
Differing  assumptions  of  the  various  calculations  are  as  indicated  in  the  figure  and  discussed  in  the 
text. 
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by  Figs.  1  and  2.  A  frequently  used  approximation  for  line  transfer  is  to  multiply  the  spontaneous 
decay  rate  by  the  escape  probability,  to  obtain  an  effective  decay  rate  which  compensates  for 
re-absoipdon.  The  results  of  Ref.  25  may  be  used  to  obtain  a  planar  average  escape  probability 
Pe  of  8.9  x  10~3  for  this  case.  Including  the  effects  of  collisional  depletion  of  the  upper  level,  the 
total  effective  decay  rate  is  ~  Aui(Pe  +  e)=  1.05  x  1012  sec-1,  implying  a  lifetime  of  0.95  psec. 
This  value  is  veiy  close  to  the  average  e-fold  time  observed  in  Figs.  1  and  2.  Transient  delays  in 
re-excitation  result  in  an  initially  more  rapid  population  falloff,  but  this  phenomenon  damps  out 
over  7  ps,  resulting  in  good  correspondence  between  the  analytically  estimated  and  numerically 
computed  upper  level  decay  times. 

The  time  profile  of  escape  of  the  He-a  line  photons  is  illustrated  in  Fig.  3  where  numerical 
results  are  plotted  along  with  the  analytic  solution  of  eqs.  (8)  and  (9).  As  indicated  in  that  figure, 
one  of  the  numerical  calculations  ignored  the  time-dependent  term  in  eq.  (1),  therefore  implicitly 
neglecting  the  photon  travel  time.  The  analytic  solution  also  neglects  this  effect.  The  numerical 
solutions  which  do  include  full  time  dependence  were  obtained  for  both  2  and  4  Gaussian  rays, 
and,  as  evident  in  Fig.  3,  die  agreement  is  excellent.  All  results  presented  in  succeeding  sections 
were  therefore  obtained  using  the  more  economical  2  ray  approximation. 

The  asymptotic  or  ultimate  escape  probabability  =  P«(  jff ) is  0.23  for  the  escape  and 
quenching  probabilities  quoted  above.  Most  of  the  emitted  He-a  photons  do  not  escape,  since  the 
quenching  probability  (e  =  0.029)  exceeds  that  of  escape  (Pe  =  8.9  x  10-3)  by  a  factor  of  3.  The 
frilly  time  dependent  numerical  solutions  depart  from  the  analytic  result  in  several  details.  First, 
the  photons  do  not  emerge  immediately  since  170  fsec  is  required  to  traverse  the  50  fim  from  the 
center  to  the  outer  boundary  of  the  plasma.  Also,  about  28%  of  the  initially  created  line  photons 
escape  rather  than  the  23%  g  ven  by  the  analytic  solution.  This  is  not  particularly  surprising 
since  the  fixed  escape  probability  is  computed  from  the  plasma  midplane  but  the  photons  diffuse 
spatially  toward  the  boundary  and  therefore  should  not  strictly  obey  a  single  escape  probability. 
According  to  the  fully  time  dependent  numerical  calculation,  1 .3  psec  is  required  for  50%  of  the 
photons  to  escape.  Achieving  a  subpicosecond  x-ray  pulse  in  A1  He-a  may  be  difficult  for  a  plasma 
like  the  present  one  where  the  line  scatters  many  times  prior  to  escape.  Note  that  recombination 
of  the  heliumlike  species  occurs  on  a  timescale  of  >  10  psec  [5],  even  for  electron  densities  of  ~ 
1022  cm-3.  Therefore  it  will  not  eliminate  the  source  of  the  multiple  scattering  during  the  x-ray 
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.  Escaped  photon  fraction  for  the  A1  He-a  line  is  illustrated  for  the  conditions  of  Figs.  1  and 
function  of  time.  Numerical  calculations  using  4  Gaussian  rays  assuming  both  the  con-ect 
finite  light  velocity  are  shown,  as  well  as  2  ray  and  analytic  approximations  discussed  in  the 


pulse.  For  near  solid  density  plasmas,  three-body  recombination  is  fast  enough  to  short-circuit 
this  mechanism.  More  general  conditions  and  assumptions  are  now  examined,  and  more  favorable 
plasma  properties  conducive  to  shorter  pulses  are  considered. 

B.  Multilevel  aluminum  plasma 

In  this  subsection,  the  simplifying  approximation  of  a  parameterized  two-level-atom  is 
dropped  and  the  escape  and  diffusion  of  the  He-a  line  is  examined  with  a  multistage,  multilevel 
atomic  model  for  aluminum  ions.  Due  to  the  expense  of  such  fully  time-dependent  calculations, 
it  was  not  feasible  to  do  a  truly  comprehensive  set  of  runs  closely  covering  all  of  the  expected 
conditions  where  such  effects  might  be  important.  Nevertheless,  with  the  aid  of  the  analytic 
model  discussed  above  and  the  atomic  number  scaling  considerations  presented  below,  essential 
conclusions  can  be  drawn  for  much  of  the  important  physics  as  applied  to  frequently  achieved 
conditions  in  laboratory  experiments. 

The  aluminum  atomic  model  emphasizes  the  heliumlike  stage  where  excited  levels  are  carried 
ranging  through  n=7.  Triplet  and  singlet  sublevels  are  modeled  for  n=2  and  n=3,  but  n=4  through 
n=7  are  treated  as  Rydberg  levels  in  which  the  triplets  and  singlets  are  assumed  to  be  populated 
statistically.  A  total  of  33  ground  and  excited  states  are  included  in  the  neutral  through  bare  nuclear 
stages.  The  heliumlike  rates  affecting  the  n=2  sublevels  are  of  the  most  importance  for  the  present 
study  and  are  in  good  agreement  with  those  of  Zhang  and  Sampson  [33].  For  the  methods  used  for 
the  atomic  rate  calculations  the  reader  may  consult  Ref.  20  and  references  therein. 

The  profile  of  the  A1  He-a  line  is  assumed  to  be  described  by  the  Voigt  function  with  the 
broadening  parameter  obtained  by  summing  the  rates  out  of  the  upper  and  lower  levels.  This 
line  is  one  of  the  last  transitions  in  the  heliumlike  manifold  to  be  affected  by  Stark  broadening  as 
density  increases.  Detailed  calculations  of  the  Static  profile  of  He-a  have  been  presented  in  Ref. 
34,  for  neon  and  argon,  for  typical  temperatures  at  which  these  ions  exist.  For  neon  (Z=10),  the 
Stark  profile  is  less  than  10%  of  the  Dojjpler  width  for  electron  densities  <  1022  cm-3.  In  the 
case  of  argon  (Z=18),  electron  density  of  at  least  2  x  102S  is  needed  for  the  Stark  width  to  exceed 
10%  of  the  Doppler  width.  According  to  Ref.  35,  the  Stark  width  in  eV  scales  as  N,/3Z-1 .  The 
temperature  at  which  the  heliumlike  stage  exists  is  approximately  proportional  to  Z3  [33],  the  ionic 
mass  scales  as  Z,  and  line  energy  as  Z2,  leading  to  a  line  Doppler  width  scaling  of  Z3.  Therefore, 
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the  electron  density  at  which  the  line  Stark-to-Doppler  width  ratio  is  a  constant  0.1  would  scale 
as  Z6 .  The  detailed  calculations  quoted  above  are  in  reasonable  accord  with  this  approximation  at 
Z5*1 .  It  is  concluded  that  for  aluminum,  the  Stark  width  of  the  He-a  line  would  not  exceed  10%  of 
the  Doppler  width  except  for  electron  densities  >  4  x  1022  cm-3.  At  such  densities,  the  collisional 
quenching  parameter  exceeds  0.5,  and  the  line  photons  can  no  longer  undergo  multiple  scattering. 
Therefore,  for  present  purposes,  the  Voigt  profile  approximation  is  reasonable. 

Results  of  the  detailed  configuration  multilevel  calculations  for  aluminum  are  shown  in  Figs. 
4-7.  In  Figs.  4-6,  the  total  ion  density  is  assumed  to  be  3.0  x  101#,9.1x  101#,and9.1  x  lO^cm-3, 
respectively.  A  calculation  at  solid  density  (6  x  10 22  cm-3)  was  also  done  and  is  discussed  below. 
The  plasma  thickness  is  100  y m  in  Figs.  4  and  5,  but  10  ym  in  Fig.  6.  In  each  case  the  plasma 
temperature  is  200  eV,  but  there  is  no  reason  to  expect  qualitative  differences  for  any  temperature 
at  which  the  heliumlike  stage  dominates  (i.e.,  150-400  eV  for  Al).  As  is  evident  from  die  “t=0” 
curves  in  these  figures,  the  initial  conditions  ate  essentially  the  same  as  for  the  two-level-atom 
case.  Initially,  the  innermost  10%  of  the  planar  plasma  has  the  1S2P1  Pi  level  excited  to  a  fractional 
population  of  0.01.  The  spread  of  excitation  due  to  the  diffusion  of  the  He-a  line  photons  is 
apparent.  Modification  of  the  escape  probability  due  to  the  Doppler  shift  of  plasma  bulk  motion 
was  not  taken  into  account;  such  effects  are  discussed  in  the  next  subsection. 

As  the  electron  density  increases  in  Figs.  4-6,  die  collisional  quenching  probability  per 
scattering  e  increases  from  0.01  (Fig.  4)  to  0.029  (Fig.  5)  to  0.23  (Fig.  6).  The  single  flight 
photon  escape  probabilities  Pe  for  these  three  cases  are  1.6  x  10~2,  8.9  x  10~3,  and  1.5  x  10~2, 
respectively.  Even  though  the  column  density  of  plasma  represented  in  Figs.  5  and  6  is  the  same, 
the  escape  probability  increases  at  the  higher  density  due  to  the  enhancement  of  the  optically 
thin  line  wings  in  the  assumed  Voigt  profile.  However,  the  ultimate  escape  probability  Pu  that  a 
photon  leaves  the  plasma  after  one  or  multiple  scatterings,  drops  monotonically  with  increasing 
density  because  of  the  dominant  influence  of  collisional  quenching.  The  expression  given  in  Sec. 
II  B  yields  P„  =  0.61,  0.23,  and  0.05  for  the  conditions  of  Figs.  4-6,  respectively.  The  differing 
collisional  destruction  of  the  He-a  photons  is  clearly  evident  in  Figs.  4-6.  The  spatial  spread 
of  excitation  with  time  is  much  more  pronounced  in  the  lower  density  cases  of  Figs.  4  and  5. 
whereas  at  the  ion  density  of  9.1  x  1020  cm-3  represented  by  Fig.  6,  it  is  very  difficult  for  a  typical 
He-a  photon  to  run  the  gauntlet  of  collisional  destruction  and  reach  the  plasma  boundary.  The 
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FIG.  4.  A1  XII  ls2p>Pl  population  fraction  is  shown  as  a  function  of  space  and  time  as  obtained 
from  a  multilevel  fully  time-dependent  coupled  rate  and  radiative  transfer  calculation  as  described 
in  the  text.  The  planar,  100  fim  thick  plasma  has  an  assumed  temperature  of  200  eV  and  total  ion 
density  of  3.0  x  101'  cm-*. 
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FIG.  5.  As  in  Fig.  4,  except  that  the  total  ion  density  is  9.1  x  10'*  cm-5,  yielding  an  electron 
density  of  10J1  cm-*. 
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conditions  of  Fig.  5  are  nearly  identical  to  those  for  which  the  parameterized  two-level-atom  case 
(Figs.  1-3)  was  studied  in  the  preceding  section.  The  excitation  profiles  as  a  function  of  time  for 
the  multilevel  and  two-level-atom  cases  may  be  compared  in  Figs.  1  and  5  at  0.5  and  2.0  ps.  They 
show  very  good  qualitative  agreement;  exact  correspondence  is  not  expected  since  quenching  of 
the  line  is  affected  by  multilevel  kinetics. 

Fig.  7  quantifies  the  time  dependence  of  the  escape  of  the  He-a  line  radiation  for  these  three 
cases  as  calculated  with  the  detailed  multilevel  model.  At  t=5  ps,  nearly  all  the  photons  have 
escaped  and  the  values  of  each  of  the  curves  of  Fig.  7  may  be  taken  as  a  good  estimate  of  Pu  •  They 
indicate  that  the  analytic  model  predictions  quoted  above  are  quite  respectable  estimates:  0.61, 
0.23,  and  0.05  vs.  the  computed  values  of  0.64, 0.27,  and  0.07.  Of  course,  the  exact  escape  time  is 
partly  a  matter  of  definition,  but  half  the  photons  which  ultimately  will  escape  have  escaped  within 
~  0.3  ps  for  the  highest  density  case,  ~  1.6  ps  for  the  medium  density  of  9.1  x  1019  cm-3,  and 
within  ~  1 .8  ps  at  the  lowest  density  considered.  At  high  density,  the  photons  are  quenched  if  they 
do  not  escape  the  plasma  on  first  flight  or  after  a  few  scatterings,  hence  the  short  pulsewidth.  At 
the  solid  ion  density  of  6.0  x  1022  cm-3,  the  quenching  probability  for  A1  is  0.95,  and  a  detailed 
calculation  confirmed  the  obvious  expectation  that  the  photon  either  escapes  after  initial  collisional 
excitation  of  the  ls2p1P1  level,  or  is  destroyed.  No  spatial  diffusion  occurs  in  this  situation  and 
no  lower  limit  on  the  pulsewidth  is  set  by  these  effects.  The  higher  level  resonance  lines  of  the 
heliumlike  ion  are  of  little  significance  for  the  present  considerations.  They  are  nearly  always 
much  weaker  spectroscopically  than  the  principal  n=2  resonance  line.  Also,  the  higher  lines  are 
much  more  susceptible  to  collisional  quenching  because  the  radiative  decay  rates  are  lower  and 
the  collisional  rates  from  the  upper  levels  are  much  higher  than  those  of  He-a. 

C.  Generalized  plasma  conditions 

1.  Effect  of  bulk  motion 

It  is  well  known  that  the  frequency  Doppler  shift  of  the  line  profile  due  to  expansion  or 
compression  of  a  medium  can  significantly  enhance  the  escape  probability.  For  planar  geometry 
(to  which  the  present  work  is  confined),  Pe  may  be  calculated  according  to  the  method  of  Sobolev 
[28]  when  the  velocity  gradient  is  constant.  Ref.  29  extends  Sobolev’s  work  to  cylindrical  media. 
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FIG.  7.  Escaped  photon  fraction  is  plotted  as  a  function  of  time  for  the  A1  He-a  line.  Total  A1  ion 
densities  and  other  plasma  conditions  are  as  assumed  in  Fi»«  4-fi 


The  frequent  occurence  of  regions  of  nearly  constant  velocity  gradient  in  numerical  and  analytic 
solutions  of  plasma  hydrodynamic  behavior  has  fostered  widespread  use  of  these  methods  to 
obtain  line  escape  probabilities.  A  full  numerical  treatment  would  require  solution  of  the  radiative 
transfer  equation  in  the  comoving  frame.  However,  the  analytic  model  embodied  in  eqs  (5)  and 
(6)  needs  no  modification,  only  the  substitution  of  the  Sobolev  escape  probability  for  the  static  one. 
The  effect  of  motional  Doppler  shifts  may  already  have  been  detected  in  the  observed  broadening 
of  heliumlike  A1  lines  [36]. 

The  magnitude  of  the  effect  may  be  estimated  using  the  analytic  solution  for  isothermal 
plasma  expansion  discussed  in  Ref.  37.  The  ion  acoustic  velocity  fox  Ai  is  ~  107  cm  sec-1  at  200 
eV.  A  1  nsec  prepulse  would  allow  plasma  expansion  to  ~  100  /xm,  implying  a  velocity  gradient 
of  10°  sec-1.  Two  examples  of  100  /xm  planar  Al  plasmas  were  discussed  above  in  the  static 
limit.  For  the  low  density  case  of  3.0  x  101®  ions  cm-3 ,  the  motionally  induced  escape  probability, 
calculated  from  eq.  (40)  of  Ref.  28,  is  8.0  x  10~3,  compared  to  1.6  x  10~2  for  the  static  limit.  At 
the  higher  density  of  9.1  x  101®  cm-3,  the  Sobolev  escape  probability  is  2.6  x  10-3,  as  opposed  to 
the  static  value  of  8.9  x  10~3.  At  this  assumed  velocity  gradient,  therefore,  the  escape  probability 
is  at  most  only  slightly  enhanced  by  motion  and  little  effect  is  anticipated. 

However,  higher  velocity  gradients  than  10®  sec-1  are  anticipated  in  some  instances, 
especially  for  kilovolt  temperature  plasmas,  as  seen  in  the  simulations  of  Ref.  36  where  values 
of  2  x  1010  sec-1  are  calculated  within  1  psec  of  the  main  pulse.  Applying  this  much  higher 
velocity  gradient  to  the  two  plasmas  discussed  above  increases  Pz  to  0.16  and  0.053,  at  densities 
of  3.0  and  9.1  x  101"  cm-3,  respectively.  According  to  eq.  (10),  the  time  for  half  of  the  He-a 
photons  to  escape  would  be  reduced  by  factors  of  5  and  2  for  the  two  cases,  respectively.  It  is 
therefore  possible  that  mo/^onal  broadening  effects  may  significantly  shorten  the  pulsewidth  of  a 
strong  x-ray  line,  depending  on  the  details  of  the  particular  plasma’s  hydrodynamics. 

2.  Atomic  number  scaling 

The  potential  for  broadening  of  the  pulsewidth  of  a  line  is  greatest  when  both  the  escape 
and  quenching  probabilities  are  small,  allowing  many  scatterings  within  the  plasma  In  this  limit, 
according  to  eq.  (10),  the  escape  time  scales  as 

fe  -  [-4„,(Pe  +0]~!  (11) 
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Given  current  experimental  capabilities,  the  atomic  number  region  8  <  Z  <  20  is  of  most  interest 
for  production  of  significant  K  shell  emission.  The  spontaneous  decay  rate  for  He-a  has  been 
calculated  and  tabulated  as  a  function  of  Z  in  Ref.  1;  the  relevant  collisional  rate  coefficients  for 
the  n=2  levels  may  be  found  in  Ref.  33.  The  Einstein  coefficient  Aui  scales  approximately  as 
Z4-2;  the  collisional  rate  coefficients  as  Z-3,2,  leading  to  a  very  sharp  dependence:  c  ~  Z-7*4.  The 
collisional  quenching  probability  for  Si  He-a,  for  instance,  is  only  about  half  that  for  an  identical 
A1  plasma,  with  a  difference  of  only  one  atomic  number.  The  scaling  of  escape  probability  Pt 
is  considerably  weaker.  According  to  eq.  (11),  the  faster  decay  rate  for  higher  atomic  number 
elements  partially  compensates  the  much  lower  collisional  quenching,  but  (neglecting  the  weak 
scaling  of  Pe),  the  escape  time  can  increase  with  atomic  number  as  fast  as  Z3'2.  Stated  another  way, 
plasma  density  must  be  increased  along  with  atomic  number  to  maintain  the  ratio  of  collisional  to 
radiative  processes  and  thereby  prevent  additional  broadening  of  the  line  photon  emission  pulse. 

IV.  CONCLUDING  REMARKS 

In  plasmas  containing  mostly  K  shell  ionization  stages,  the  strongest  resonance  line  (He-a) 
is  usually  a  significant  if  not  dominant  component  of  the  x-ray  emission.  This  is  true  of  plasmas 
lasting  several  nsec,  such  as  the  dense  Z  pinch,  as  well  as  of  plasmas  produced  by  ultrashort  pulse 
lasers  which  may  persist  for  psec  timescales  or  shorter.  Photons  emitted  in  this  line  as  well  as 
other  strong  resonance  lines  can  be  subject  to  many  absorptions  and  re-emissions  in  optically  thicx 
plasmas.  In  general  this  time-broadens  the  line  component  of  the  x-rays. 

This  paper  has  presented  detailed  numerical  calculations  supplemented  and  interpreted  by 
an  analytic  model  which  demonstrate  that  for  laboratory  plasmas  of  moderate  density  and  mm 
or  smaller  size,  the  escape  time  is  a  few  psec  or  less.  These  models  could  also  be  employed  in 
analyzing  astrophysical  scale  plasmas  where  much  longer  characteristic  times  are  expected.  One 
current  goal  of  experiments  employing  high  power  subpicosecond  lasers  is  to  produce  the  shortest 
possible  x-ray  pulse.  If  the  entire  laser-produced  plasma  is  of  solid  density,  collisional  quenching 
of  the  line  photons  assures  that  the  scattering  mechanism  plays  no  role  in  widening  the  x-ray  pulse. 
However,  to  date  the  most  efficient  x-ray  generation  has  occured  when  a  preformed  plasma  created 
by  a  laser  prepulse  is  subsequently  irradiated  by  the  main  high  intensity  pulse.  The  persistence 
of  the  heliumlike  ground  state  raises  the  possibility  that  a  photon -scattering  surrounding  plasma. 
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even  though  much  cooler  than  the  principal  x-ray  generating  region,  may  place  a  lower  limit  on  the 
width  of  the  x-ray  pulse.  Perhaps  the  use  of  mostly  continuum  emission  or  lines  of  other  ionization 
stages  which  may  be  more  subject  to  collisional  quenching  would  ameliorate  this  effect.  The 
enhancement  of  the  line  escape  probability  due  to  the  motional  Doppler  shift  may  further  reduce  its 
importance  in  some  instances.  The  time-dependent  transfer  of  radiation  in  strong  resonance  lines 
is  one  of  several  physical  processes  which  determine  the  shape  and  intensity  of  the  x-ray  emission 
from  the  new  class  of  ultrashort  lived  laser-produced  plasmas. 
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Abstract 


If,  through  stimulated  Brillouin  scattering  or  through  parametric  laser  absorption  near  the 
critical  surface,  a  significant  amount  of  ion-acoustic  microturbulence  is  generated  in  laser-produced 
plasmas,  then  a  10-  to  100-fold  reduction  in  the  intensity  threshold  for  generating  non-Maxwellian 
distributions  of  the  kind  described  by  Langdon  (A  B.  Langdon,  Phys.  Rev.  Lett  44, 575  (1980)) 
could  occur.  Because  these  distributions  have  depleted  tails,  they  would  alter  the  ionization 
dynamics  of  these  plasmas  and  significantly  affect  the  interpretation  of  x-ray  laser  experiments,  in 
particular.  In  this  paper,  we  will  discuss  the  potential  effects  of  laser-produced  non-Maxwellians 
on  the  gain  dynamics  of  selenium  neon-like  x-ray  lasers.  We  utilize  an  atomic  model  with  a 
full  set  of  n=3  neon-like  multiplet  states,  and  compare  non-Maxwellian  J=0-1  and  J=2-l  gain 
calculations  to  corresponding  Maxwellian  calculations.  It  is  found  that  die  generation  of  a  strong 
ion-acoustic  turbulence  at  laser  intensities  in  excess  of  1015  W/cm2  could,  in  principle,  help 
to  promote  increases  in  the  gain  of  the  3s-3p  transitions  in  neon-like  selenium  and  to  promote 
significant  improvements  in  the  performance  of  these  lasers. 
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I  Introduction 

In  laser-produced  plasmas,  there  are  many  problems  associated  with  the  absorption  of 
laser  light  in  the  plasma.  In  this  case,  the  absorption  of  electromagnetic  energy  by  inverse 
Bremsstrahlung  can  be  strongly  perturbed  by  the  competing  processes  of  stimulated  Raman  or 
stimulated  Brillouin  scattering.  These  scatterings  can  generate  a  strong  microturbulence  through 
the  nonlinear  generation  of  plasma  waves  and  ion-acoustic  waves.1  When  laser  absorption  occurs 
at  the  critical  surface,  moreover,  parametric  decay  instabilities  are  excited  that  generate  a  strong 
Langmuir  turbulence  that  decays  into  ion-acoustic  turbulence.  Strong  temperature  gradients 
generated  by  the  absorption  of  laser  light  in  sub-critical  x-ray  laser  plasmas  can  initiate  electrical 
currents  of  sufficient  strength  to  generate  ion-acoustic  turbulence.  Unfortunately,  it  is  difficult  to 
calculate  the  absorption  coefficient  of  laser  light  in  microturbulent  plasmas  unless  die  dynamical 
growth  and  saturation  of  the  microturbulence  energy  spectrum  can  be  calculated  throughout  die 
plasma. 

Microturbulenoe  also  affects  the  heat  transport  and  viscosity  properties  of  a  collisional  plasma 
as  well  as  the  stupe  of  die  isotropic  electron  distribution  function.2'3  These  microturbulence 
properties  will  manifest  themselves  in  a  variety  of  ways  in  laser-produced  plasmas,  and  die  full 
combination  of  their  effects  should  be  used  to  diagnose  the  presence  of  microturbulence.  Moreover, 
microturbulence  may  need  to  be  factored  into  any  analysis  of  laser-produced  plasma  dynamics 
before  other  effects  such  as  nonlocal  heat  transport  can  be  invoked4  to  fully  explain,  for  example, 
the  inhibited  heat  flow  known  to  be  present  in  these  plasmas.5 

In  this  paper,  we  assume  that  a  fully  developed  microturbulence  has  been  generated  in  a  laser- 
produced  selenium  plasma.  Depending  on  the  charge  states  of  the  plasma,  the  microturbulence 
will  change  the  plasma’s  laser  absorption  as  well  as  its  transport  properties.  In  this  paper,  we 
will  calculate  selenium’s  charge  state  in  Collisional-Radiative  Equilibrium  (CRE)  utilizing  the 
atomic  model  described  in  Ref.  (6)  in  order  to  illustrate  the  effect  that  the  non-Maxwellian 
electron  distributions  generated  by  the  laser  absorption  process  can  have  on  the  neon-like  gain 
coefficients.  When  ion-acoustic  turbulence  is  generated  as  a  result  of  the  laser-plasma  interaction, 
this  microturbulences  makes  it  easier  to  generate  these  non-Maxwellian  distributions. 

We  begin  by  initially  investigating  the  coupling  of  ion-acoustic  turbulence  to  the  electron 
kinetic  equation.  We  then  seek  to  determine  the  plasma  conditions  and  ion-acoustic  fluctuation 
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levels  needed  to  significantly  enhance  a  weakly  turbulent  selenium  plasma’s  laser  absorptivity 
given  a  CRE  charge  state  of  the  plasma  as  a  function  of  electron  temperature  and  ion 
density.  Microturbulence  promotes  the  production  of  self-similar  non-Maxwellian  electron 
velocity  distributions  of  the  kind  described  by  Dum2  and  Langdon,7  i.e.,  Langdon’s  criterion  for 
non-Maxwellian  production  is  modified  by  the  presence  of  a  fully  developed  microturbulence.  The 
laser  intensities  that  are  needed  to  drive  selenium  into  non-Maxwellian  states  are  then  determined 
from  this  criteria.  These  results  may  have  special  significance  to  x-ray  laser  research.  Neon¬ 
like  selenium  is  a  well  studied  x-ray  laser  medium,8  and  a  number  of  discrepancies  regarding 
its  observed,  versus  its  calculated,  ionization  behavior  remain  to  be  explained.®  Some  of 
this  discrepancy  might  be  explained  by  the  presence  of  microturbulence  and  its  promotion  of 
non-Maxwellian  ionization  behavior. 

This  paper  is  structured  as  follows.  Section  II  contains  a  new  derivation,  more  general 
than  one  presented  by  Dum,3  of  the  way  ion-acoustic  turbulence  can  be  included  in  the  kinetic 
theory  of  a  collisional  plasma.  Separate  contributions  of  this  microturbulence  to  the  different 
terms  of  the  expansion  of  the  distribution  function  in  spherical  harmonics  are  derived.  These 
contributions  are,  to  good  approximation,  identical  in  form  to  the  ones  made  by  election-ion 
collisions.  Thus,  Langdon’s  criterion  can  easily  be  rederived  for  the  laser  intensity  needed  to 
drive  the  elections  into  a  non-Maxwellian  state  due  to  rapid  heating  in  a  microturbulent  plasma.7 
Using  a  phenomenological  model  for  the  fluctuation  level  scaling  achievable  in  weakly  turbulent 
plasmas,10  we  then  calculate  the  selenium  plasma  conditions  under  which  microtuibulence  would 
be  expected  to  produce  non-Maxwellian  electron  distributions  as  a  function  of  selenium  plasma 
conditions.  In  section  HI,  the  effects  of  these  non-Maxwellians  on  gain  calculations  in  neon-like 
selenium  is  presented  and  discussed.  Finally,  some  implications  of  this  work  for  experimentally 
diagnosing  plasma  microturbulence  and  for  improving  the  performance  of  neon-like  selenium 
lasers  are  discussed  in  section  IV. 

II  Microturbulence  Analysis 

We  begin  our  analysis  by  writing  the  electron  distribution  function,  /e,  as  the  sum  of  an 
average  (in  space  and  time),  {/),  and  of  a  rapidly  fluctuating  distribution,  6f:  fe  =  {/)  +  £/. 11  We 
work  in  the  coordinate  system  of  the  fluid  flow,  Fourier  transforming  the  spatial  coordinates  from 
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x  to  k,  and  utilizing  standard  quasilinear  theory  to  solve  for  if  by  ignoring  particle  accelerations 
and  plasma  fluid  flow:11 


«/(k,v) 


“ _ 53H _ v.</> 

m  u;(k)  —  k  •  v  +  *7(k)  ’ 


(1) 


where  e  is  the  electronic  charge,  a  positive  number,  and  m  is  the  electron  mass.  The  time 
dependence  of  the  electromagnetic  fluctuations  is  given  by  exp[t(o>(k)  +  iy(k))tj.  This  expression 
for  if  must  be  substituted  into  the  equation  for  {/),  which  contains  a  collision  term,  C((f)),  that 
is  the  sum  of  electron-electron,  electron-ion,  and  electron-microturbulence  interactions: 


C  =  C“  +  Cei  +  CAN , 


(2) 


where 

CA"  s  =  -<«E(r,/).V.*/(r,v,«)).  (3) 

m 

Note,  because  we  have  specialized  to  ion-acoustic  microturbulence,  the  more  general  notation 
CAN  will  be  replaced  by  Cia.  When  Eq.  (1)  is  substituted  into  Eq.  (3),  a  result  that  can  be  written 
in  terms  of  a  microturbulence  diffusion  tensor,  is  obtained: 
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where 
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(5) 

(6) 


W(k)  is  the  energy  spectrum  of  the  ion-acoustic  fluctuations,  i.e.,  (|f.E(r,<)|2/87r)  =  J  dsfcW(k), 
and  ki  =  k,/Aj.  We  are  interested  in  evaluating  1 for  those  physical  situations  in  which  the 
turbulence  is  fully  developed  and  the  plasma  is  marginally  stable.12  In  these  situations,  the 
limit  7  — ♦  0  can  be  taken  so  that  7 /((a;  —  k  •  v)2  -f  y2)  — >  7r£(u;  —  k  •  v).  We  further  simplify 
the  problem  by  assuming  that  the  microturbulence  is  isotropic.  (The  more  complete  theory  of 
non-isotropic  microturbulence  is  reviewed  by  Bychenkov,  Silin,  and  Uryupin.13) 

The  dispersion  relation  for  ion-acoustic  waves  is  u>(k)  =  ±\/k2c2J{  1  +  k7 X2D)  =  ±kct, 
where  c,  =  J~ZkBTt/mi  is  the  ion  sound  speed,  and  A d  =  vth/(x/2we)  is  the  electron  Debye 
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length.  vth  =  y/2 kaTe/m  is  the  average  electron  thermal  velocity,  we  =  ^/4irnee2 /m  is  the 
electron  plasma  frequency,  Z  is  the  average  charge  state  of  an  ion,  and  m»  is  the  ion  mass.  The 
diffusion  tensor  can  now  be  evaluated  by  letting  the  spherical  coordinate  unit  vectors,  ug,  u^,  and 
uVf  in  velocity  space  define  the  directions,  kt,  kv,  and  kz,  in  k-space  respectively.  The  k-space 
integration  in  Eq.  (S)  can  be  carried  out  in  the  spherical  k-space  coordinates,  (k,  O' ,  </>'),  where 
k  *  v  =  kv  cos  6'  =  kvfi.  Using  the  reality  condition,  u>(-k)  =  -u>(k),  we  rewrite  Eq.  (5)  as 


g_  2-2  roo  f  *0 

Vij  =  -—j-  j  k2dkW(k)J  J  dfik{kj6(kvfi  +  u(k)) 

+  J  dfikikj£(kvfi  —  w(A:))|,  (7) 

where  u>(k)  =  kc,  for  ion-acoustic  waves.  The  integrations  produce  the  following  diagonal  diadic 
form: 


V  =  V0(ugug  +  u^u^)  4-  (Vo  +  Vi)uvuv. 

(8) 

For  ion-acoustic  turbulence. 

^  4ir2e2Hi 

Vo  —  , 

mrv 

-  jj  j"7+(»  -  c.), 

(9) 

47r2e2Hi 

X?i  = - - — - 

mrv 

j1' 

(10) 

In  these  expressions. 

fOO 

III 

H 

CJ 

/  dkknW(k), 

'o 

(11) 

V+(x)  is  the  Heaviside  function,  which  equals  one  when  *  >  0  and  is  zero  otherwise.  The  spectrum 
integration,  fli  defines  the  strength  of  the  ion-acoustic  microturbulence.  Note  that  ion-acoustic 
turbulence  affects  essentially  the  entire  distribution  function,  because  c,  <C  vtfl,  and  one  can  take 


V0  =  -Pi  = 


47rJe2Hi 


to  good  approximation  in  order  to  describe  these  interactions. 

We  next  make  the  usual  approximate  expansion  of  (/)  in  spherical  harmonics14  in  the 
coordinate  system  moving  with  the  fluid  velocity,  v/: 


(/)  =  /o  +  fi  • 


(13) 
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The  collision  terms  are  similarly  expanded:  C  =  C0  +  Cx  •  v/v.  The  equations  of  motion  for  /0 
and  fi  are  then  (see  Ref.  (15),  pages  1 14,1 15) 


(a,  +  V/  •  v)/, .  (v  •  V/):5£  +  •  f,  +  (v2ar  •  fl)  =  <w.) 


and 

( &t  4-  v/  •  V)fj  -  (Vv/)  •  fi  -  (Vv/  +  (Vv/)T  +  I2V  •  v/)  •  ~ f ) 

5  ov  \  v  J 

-"c  X  fi  -  Ci(f,)  =  — «V/o  -  a'^,  (15) 

where 

a'  “ --(E  +  —  x  B)  =  E'.  (16) 

m  c  m  v  ; 

E'  is  the  effective  electric  field  seen  in  the  local  frame  of  the  moving  fluid,  and  we  =  eB /me. 

Both  Co  and  Ci  contain  contributions  from  C“,  which  we  denote  (7j°  and  Cj"  respectively.  From 

Eqs.  (4)  and  (6),  one  finds  that 


Cj-  =  ~a.{v2(v0+vI)a.f0) 


(17) 


and 

C?  =  ^d9(v2(V0  +  Z>i)0,fi)  -  ~V0fi.  (18) 

(Similarly,  by  going  one  more  term  in  the  spherical  harmonic  expansion,  one  obtains 

cj*  =  v~2dv  [v2(x>o  +  Vi)dvr2 ]  -  6v-2v0t2.) 

Eq.  (12)  implies  that  —  0  and 


(19) 

(20) 


where  the  effective  velocity-independent  ion-acoustic  collision  frequency,  i/*“,  is 


i/*°  = 


87r2e2ni 


Thus,  like  electron-ion  scattering,  the  interaction  of  electrons  with  ion-acoustic  waves  primarily 
affects  the  momentum  transfer  to  or  from  the  electrons. 
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We  next  adopt  the  common  procedure  ofignoring  electron-ion  collisions  in  Eq.  (14)  (C"  =  0) 
and  electron-electron  collisions  in  Eq.  (15)  (C'e  =  0).  The  isotropic  electron-electron  collision 
term,  C",  is  given  by  the  Fokker-Planck  expression.16  If  the  ions  are  Maxwellian,  the  electron-ion 
momentum  transfer  in  the  ion  rest  frame  has  the  same  form  as  the  electron/ion-acoustic  term:16 

c*  = (t)1  f"  (22) 

where  the  velocity-independent  electron-ion  collision  frequency  is  defined  by 

^  (23) 

m*v*k 

and  log  Ae  is  the  Coulomb  logarithm.  Thus, 


/nr ec 

0  =  O0 


and 

C,arCf +  =  (25) 

where 

v\  =  */«»  +  */**  =  (1  +  (26) 

A  similar  expression  for  Ci  was  derived  by  Dum,3  who  carried  out  a  direct  angle  averaging  of  Eq. 
(5). 

The  effect  of  current  and  heat  flow  on  the  shape  of  the  isotropic  distribution  function,  /0,  is 
determined  by  the  two  fi  terms  in  Eq.  (14).  The  first  term,  involving  V  •  fj ,  determines  the  effects 
of  current  and  heat  flow  on  the  electron  number  and  energy  densities  respectively.  The  second 
term,  involving  af,  represents  the  effect  of  laser  heating  on  the  electrons.  To  determine  fa ,  we  look 
for  a  solution  to  Eq.  (15).  Generally,  as  a  first  approximation,  the  gradient  and  time  derivative 
terms  can  be  dropped,  leaving  one  with  the  equation, 

-wc  Xf,  +  (1  +/?H,  (^)5  f,  =  -vVfo  -  a'|£.  (27) 

This  equation  may  be  easily  solved  for  fi  in  terms  of  /o(v).  If  a  Maxwellian  distribution  is  used 
for  /o,  simple  expressions  can  then  be  found  for  the  electric  current,  je,  and  heat  flux,  qe,  from 
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which  the  thermal  conductivity,  *e,  can  he  obtained.  Both  conductivities  will  be  reduced  from 
their  classical,  non-turbulent  values  by  the  factor,  14-/?. 

To  determine  the  size  and  the  dynamics  of  these  conductivity  reductions  in  detail,  one  has  the 
difficult  task  of  calculating  the  turbulence  spectrum,  W.  However,  since  only  one  moment,  17 1 , 
of  W  is  needed  to  determine  the  size  of  fi  for  ion-acoustic  turbulence,  it  is  possible  to  estimate 
how  strongly  a  given  level  of  electrostatic  microturbulence  will  affect  these  conductivities  by 
establishing  a  bound  on  fti .  We  begin  with  the  definition  of  /?.  From  Eqs.  (21)  and  (26), 


where  £<&  is  die  electron  thermal  energy  density: 

Etk  =  (3/2  )ri'kBTe.  (29) 

The  size  of  Hi  can  be  estimated  by  factoring  an  average  ion-acoustic  wavelength,  l/(Ar),  from  its 
integrand;  then,  17i  =  feEtk/ityt  where  /e  is  the  ratio  of  the  electrostatic  fluctuation  energy  to 
the  thermal  energy,  ({|625(r,  t)\2/Sir)/Eth)-  Thus, 


If,  by  virtue  of  the  ion-acou  .tic  dispersion  relation,  one  assumes  that  (k)  <  1/A d,  then 

P  >  Pmi*  =  VW,/E.  (31) 

For  laser-produced  plasmas,  the  product,  a >ere,  can  range  in  value  from  >  1  to  >  104;  thus, 
for  large  values  of  this  ratio,  ion-acoustic  energy  fluctuation  levels,  fs,  of  1%  can  easily  lead 
to  a  turbulence  contribution  to  the  thermal  conductivity  that  is  10  times  larger  than  that  due  to 
two-particle  collisions.  This  turbulence  level  corresponds  to  a  heat  flux  limiter  of  0. 1 ,  a  value  often 
used  to  describe  the  dynamics  of  laser-produced  plasmas.5 

Investigating  the  non-turbulent  laser  heating  of  electrons,  Langdon7  derived  a  criterion  that 
predicted  the  intensity  at  which  inverse  Bremsstrahlung  heating  v,  r»uld  dominate  over  electron- 
electron  collisions.  He  found  that,  when  this  criterion  was  satisfied,  the  electron  distribution 
approached  a  self-similar  form  whose  high  energy  tail  was  depleted  relative  to  that  of  a  Maxwellian 
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distribution.  Earlier,  Dum2  had  found  that,  under  strong  ion-acoustic  turbulence,  the  electron 
distribution  tends  towards  the  same  self-similar  form  as  discussed  by  Langdon.  We  rederive  this 
criterion  for  a  microturbulent  plasma. 

By  substituting  the  general  solution  to  Eq.  (27)  into  Eq.  (14)  and  by  ignoring  the  vV/o  term 
in  Eq.  (27)  and  the  gradient  terms  in  Eq.  (14),  one  obtains  the  equation,  dtfo  =  H(fo  )  +  Co(fo)t 
which  is  the  equation  that  was  used  by  Langdon.  The  laser  heating  term  has  the  form, 

E(f\  47re2/  1  9  (v2  df.\ 

0  m2c  3v2  dv  \  a >2  + 1 ^  (t>t/»/t;)6  dv  )  * 

where  I  is  the  laser  intensity  and  u>  =  c/X  is  the  laser  frequency.  One  can  define  an  e-folding  time 
for  laser  heating,  rjfter,  by: 

later _  Etk _ Eth  foo\ 

"  "ti-E)u (33) 

where  klI  =  Air  J dvv2(m/2)v2H(fo)  and  k-l  is  the  laser  absorption  coefficient.  ff(/0) 
dominates  Co(fo)  when  the  laser  e-folding  time  for  heating  is  smaller  than  the  equilibration 
collision  time,  ree,  needed  to  establish  a  Maxwellian  distribution:  Tct/r^'er  >  1.  Since 
ree  =  Zt€,  one  can  write  this  criterion  as 

Wri?"r=  (34) 

3  V*h 

where  v0  is  the  peak  oscillation  velocity  of  the  electrons  in  the  high  frequency  (laser)  electric  field. 
Alternatively,  this  criterion  can  be  rewritten  as  a  condition  on  the  laser  intensity: 


W  V*K 


where,  by  definition  of  v0  and  v«&. 


3irmc3kgTe 


Icrit  — 


6irmc  ks-ie  er9 9 

2Ze2(l+0)X2  sec -cm2 


The  reduction  in  /crtt  by  the  factor,  (1  4-  /?),  corresponds  to  an  increase  in  the  laser  absorption 
coefficient,  by  the  same  amount. 

The  magnitudes  of  the  microturbulence  corrections  to  Icrit  depend  on  plasma  conditions  and 
can  be  determined  once  the  saturation  level,  /#,  is  known.  This  level  also  depends  on  plasma 
conditions  and  on  the  strength  of  the  energy  source  driving  the  turbulence.  We  investigate  this 
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dependence  by  arbitrarily  fixing  the  value  of  /e  for  one  set  of  plasma  conditions  and  by  scaling  /E 
as  described  by  Ichimaru.10  In  a  weakly  turbulent  plasma,  he  argued  that  the  turbulence  spectrum 
scales  as  g1/2,  where  g  is  the  plasma  parameter. 


'37> 

which  we  take  to  be  the  electron  plasma  parameter.  Note  that  in  strongly  turbulent  plasmas,  /e 
would  scale  as  g°,  i.e.,  it  would  be  independent  of  g.10 

To  illustrate  how  this  scaling  affects  Pmin  and  Icrit,  we  calculate  selenium’s  ionization 
state,  Z,  using  the  atomic  model  described  in  Reference  (6)  by  assuming  that  the  plasma  is  in 
CRE.  Although  laser-produced  selenium  plasmas  have  been  well  studied  recently  in  x-ray  laser 
experiments,8,9  these  plasmas  are  also  of  interest  in  z-pinch  experiments.  The  Z  values  for 
selenium  that  are  needed  to  calculate  the  quantities,  /e,  Pmin,  and  Icrit,  are  displayed  in  Fig.  (1) 
as  a  function  of  temperature  and  density.  For  the  purpose  of  this  numerical  example,  we  took 
fs  =  This  formula  produced  the  values  for  fg  that  are  shown  in  Fig.  (2)  for  the 

same  plasma  temperatures  and  densities  as  Fig.  (1).  Fig.  (2)  is  an  illustration  of  the  fact  that,  in 
weakly  turbulent  plasmas,  as  in  quiescent  plasmas,  relatively  higher  levels  of  electrostatic  energy 
fluctuations  can  be  generated  as  the  plasma  temperature  is  lowered  or  as  the  plasma  density  is 
increased.  These  fluctuation  levels  produce  the  Pmin  enhancement  factors  (see  Eq.  (31))  shown  in 
Fig-  (3),  which,  in  turn,  produce  the  contours  for  Icrit  that  are  presented  in  Figs.  (4). 

The  1^  values  in  Fig.  (4)  were  calculated  for  frequency  tripled  neodymium  (A  =  0.35 /tm) 
laser  radiation.  They  are  close  to  those  used  in  current  experiments;  thus,  as  these  experiments 
are  scaled  to  higher  photon  energies  or  higher  x-ray  laser  intensities,  the  likelihood  for  generating 
non-Maxwellian  electron  distributions  under  microturbulent  plasma  conditions  may  increase.  Note 
that  Fig.  (2)  can  be  scaled  directly  to  calculate  Pmin  in  selenium  for  any  /e  since  /?mm  is  directly 
proportional  to  /e  .  In  the  gain  calculations  that  follow,  we  will  focus  our  attention  on  the  point  at 
Te  =  1  keV  and  JV<  =  102°  ions/cm3  where  we  will  assume  that  p  =  10. 


Ill  Gain  Calculations 

When  I  »  Icrit,  the  electron  distribution  approaches  the  self-similar  form, 
u~3exp(—v5 /5us),  with  us  ~  t,  discussed  by  Dum2  and  Langdon.7  The  number  of  electrons  in 
the  tail  of  this  distribution  is  greatly  reduced,  leading  to  reductions  in  the  rates  at  which  the  plasma 
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is  ionized.  Significant  departures  of  /0  from  a  Maxwellian  distribution  also  change  the  calculation 
of  fi,  which,  in  turn,  modifies  the  electrical  and  heat  conductivities  that  are  calculated  due  to 
plasma  microturbulence  alone  (see  Refs.  (2)  and  (3)).  For  weaker  intensities,  the  distribution 
function  has  transitional  behavior  of  the  kind  calculated  by  Alaterre  Matte,  and  Lamoureux.16 
These  authors  found  that  for  laser  intensities  comparable  to  Ierit,  solutions  to  the  Fokker-Planck 
equation,  dtfo  —  H(fo)  +  Co(/o),  can  be  adequately  approximated  by  the  functions: 

=  (Cm/usm)exP  (-(•-/«*)),  (38) 


where 


Cm  = 


4flT(3/m)  ’ 


_  /  r(5/m)  \ 


'  \3T(3 

Whoa  m  =  5,  these  functions  coincide  with  the  self-similar  solutions  of  Dum  and  Langdon. 
When  m  =  2,  they  coincide  with  Maxwellian  distributions.  Thus,  for  2  <  m  <  5,  these 
functions  describe  transitional  behavior  between  rite  Maxwellians,  present  when  J  <  lent,  and 
the  self-similar  distributions,  present  when  I »  /«•*«•  The  differences  between  the  functions,  f™, 
and  a  Maxwellian  distribution  are  conveniently  represented  by  the  departure  functions,  7,  defined 

by  y(B)  =  /r(g.r.)  _  m 

iv/  —  rp  ^  fms2/r  rr  \ '  '  / 


fr~’(E,T.y 


These  functions  are  displayed  in  Fig.  (5)  for  various  m’s:  2  <  m  <  5. 

Alaterre,  et.  al.  were  able  to  fit  their  Fokker-Planck  solutions  using  the  above  functions  by 


determining  m  from  the  formula. 


m(a)  =  2  + 


1  -f-  1.66/a0-724  ’ 


where 


a  =  2 Z-%-  =  4  x  10 


_16  z I[W l cm7\\{inn\‘ 


(!+/?)• 


**  —  2  —  *  ''  rr>  r  I  t/1  r  v  / 

v2th  TeifceV] 

The  factor  1  +  /?  does  not  appear  in  their  work.  It  is  added  here  and  represents  the  natural  extension 
of  their  calculations  to  a  microturbulent  plasma.  The  transition  of  m  between  2  and  5  as  a  function 
of  laser  intensity,  when  =  10,  is  shown  in  Fig.  (6). 
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In  the  discussion  that  follows,  the  selenium  atomic  model  described  in  Ref.  (17)  is  used.  The 
model  consists  of  hydrogenic  structure  in  each  ionization  stage  of  selenium  except  the  neon-like 
stage,  which  contains  a  complete  set  of  the  26  n=3  multiplet  states  as  shown  in  Fig.  (7).  The  rate 
coefficients  for  the  L-shell  of  this  model  are  calculated  by  averaging  atomic  cross-sections  over  any 
specified  electron  distribution  function.  For  the  calculations  in  Ref.  (17),  Maxwellian  distributions 
were  used.  In  the  gain  and  ionization  abundance  calculations  that  are  presented  in  this  section, 
we  used  rates  that  were  calculated  by  integrating  over  m  =  5  distributions.  Fig.  (8)  illustrates 
the  relative  sensitivity  of  3  pairs  of  rate  coefficients  to  the  underlying  distribution  functions,  /™, 
that  are  used  to  obtain  them.  Using  Fig.  (6),  we  have  plotted  the  deviations  from  detailed  balance 
of  these  rates  as  a  function  of  I  =  IiMr  rather  than  m.  As  the  curves  in  Fig.  (5)  indicate, 
the  depletion  of  the  tail  of  the  f™  distributions  reduces  the  ionization  rate  between  the  neon-  to 
fluorine-like  ground  states.  Under  these  same  plasma  conditions,  however,  the  increases  in  the 
bulk  of  the  distributions  act  to  enhance  the  collisional  excitation  rates  to  the  J  —  0  3p  excited 
state  (denoted  CE(J  =  0)).  The  J  =  0  state  whose  rates  are  being  compared  in  this  figure  is  the 
uppermost  of  the  2p63p  multiplet  states  that  are  drawn  to  scale  together  with  the  2ps3a  states  in 
Fig.  (9).  This  figure  also  displays  the  four  J=0-1  and  J=2-l  x-ray  laser  transitions  that  have  been 
the  focus  of  research  on  neon-like  lasers. 

Figs.  (10M12)  show,  respectively,  the  effects  that  the  m  =  5  non-Maxwellian  distribution 
have  on  the  gain  coefficients  of  the  J=0-1  and  J=2-l  transitions  (C  and  D  respectively  in  Fig.  (9)) 
and  on  the  ionization  abundances  of  the  neon-  and  fluorine-like  ground  states.  The  calculations 
were  carried  out  for  purposes  of  illustration  using  CRE  solutions  to  the  atomic  rate  equations.  Note 
that  at  Tt  =  1  keV,  the  gains  of  the  C  and  D  transitions  peak  at  Ni  ~  2  x  1020  ions/cm3.  Both 
gains  are  enhanced  significantly  by  the  non-Maxwellian  distribution;  however,  the  enhancement  is 
larger  for  the  J=0-1  transition  at  168.7^4  than  it  is  for  the  J=2-l  transition  at  206.4>1.  The  reason  is 
found  in  Fig.  (12).  The  reduction  in  the  ground-to-ground  ionization  rate,  suggested  by  Fig.  (8), 
leads  to  a  significant  increase  in  the  population  of  the  neon-like  ground  state.  The  excitation  of  this 
state  is  the  dominant  mechanism  for  pumping  the  J=0,  3p  state  and  for  creating  the  3s-3p,  J=0-1 
population  inversion.  The  J=2,  3p  state,  on  the  other  hand,  is  pumped  predominantly  from  the  3d 
states  and  is  not  so  strongly  affected  by  the  shift  in  the  ionization  abundances.  Fig.  (8)  shows  that 
the  excitation  rate  to  the  J=0  state  has  also  been  increased.  Thus,  there  are  more  states  to  pump  and 
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the  pumping  rate  has  been  increased  when  the  electrons  are  distributed  in  an  m  =  5  distribution 
rather  than  anm  =  2  distribution.  The  net  effect  is  to  preferentially  populate  the  J=0  state  more 
than  the  J=2  state.  Note,  however,  that  since  the  measured  gain  coefficients  of  these  lines  in  x-ray 
laser  experiments  is  only  of  the  order  of  5,  one  can  conclude  from  Figs.  (10)  and  (11)  that  there 
is  significant  room  for  improvement  in  the  design  and  performance  of  laboratory  neon-like  x-ray 
lasers. 

IV  Summary  and  Conclusions 

The  two  major  problems  in  plasma  microturbulence  theory  are  (1)  to  determine  the 
dynamics  of  the  electromagnetic  energy  fluctuations,  inclusive  of  their  initial  growth  rates, 
their  saturation  mechanisms,  their  propagation  properties,  and  the  dependence  of  these  properties 
on  inhomogeneous  plasma  conditions,  and  (2)  to  determine  the  effects  of  microturbulence  on 
the  electron  kinetics  and  on  the  fluid  dynamics  of  the  plasma.  While  both  of  these  problems 
are  important,  only  the  latter  one  was  analyzed  in  this  paper  under  die  assumption  that  the 
microturbulence  had  been  triggered  and  had  grown  to  a  saturated  level  that  was  some  (small) 
fraction  of  the  electron  thermal  energy.  Under  this  assumption,  we  showed  how  microturbulence 
couples  to  the  the  different  terms  in  the  spherical  harmonic  expansion  of  the  electron  kinetic 
equation  for  the  case  of  ion-acoustic  turbulence.  We  then  found,  in  agreement  with  Ref.  (3), 
that  the  essential  effect  of  electron/ion-acoustic  coupling  is  to  increase  the  electron-ion  collision 
frequency.  In  so  doing,  it  alters  three  basic  properties  of  the  plasma  (by  essentially  the  same 
correction  factor):  electrical  resistivity  is  increased,  heat  flow  is  reduced,  and  the  critical  laser 
intensities  needed  to  generate  non-Maxwellian  distributions  are  reduced.  The  latter  reduction 
occurs  as  a  result  of  the  increase  in  the  laser  absorption  coefficient.  This  absorption  change  could 
lead  to  the  use  of  probe  laser  beams  to  determine  the  fluctuation  levels  of  microturbulent  plasmas 
experimentally. 

Since  microturbulence  reduces  the  threshold  for  generating  non-Maxwellian  electron 
distributions,  its  presence  would  also  manifest  itself  indirectly  through  this  distribution’s  influence 
on  the  x-ray  emissions  generated  in  the  plasma.18  These  emissions  could,  therefore,  provide 
diagnostics  of  the  microturbulence.  Because  the  laser  intensity  threshold  needed  to  produce 
non-Maxwellian  distributions  is  significantly  reduced  by  (ion-acoustic)  microturbulence,  our 
calculations  suggest  that  these  emission  effects  might  be  seen  in  laser-produced  plasma,  x-ray 
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laser  experiments.  More  importantly,  if  non-Maxwellian  distributions  are  generated  in  x-ray  laser 
experiments  conducted  with  high  atomic  number  plasmas,  then  they  could  produce  shifts  in  the 
ionization  balance  and  increases  in  pumping  rates  for  generating  population  inversions.  These 
increases  could,  in  turn,  significantly  affect  the  design  of  these  experiments  and  significantly 
improve  the  performance  of  laboratory  neon-like  lasers. 
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function  of  electron  temperature  and  ion  density 


l°9io{0mi 


Figure  3.  Contours  of  logiO(0min)  are  shown  as  a  function  of  electron  temperature  and  i 
density.  They  were  calculated  using  the  values  for  7  and  fe  shown  in  Figs.  (1)  and  (2). 
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Figure  4.  Contours  of  log 
density  using  the  Z  and  /?  vali 
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Figure  0.  T 
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Figure  7.  Atomic  structure  of  the  neon-like  ionization  stage.  All  of  the  „  =  3  mmdplet  states 
are  described  in  full  detail  within  this  atomic  model. 


DETAILED  BALANCE  OF  RATES 


XV w  / 


(oX9  —  dpA\/ 


XVW-NON  /  . 

(  OX3  —  9pA\l  OXsA\ 


lb 


Figure  8.  The  degree  to  which  rate  coefficients  faU  out  of  detailed  balance  as  a  function  of 
lnser  intensity  (or  m,  see  Fig.  (6))  is  shown.  Cl  denotes  collisional  ionization  and  recombination 
processes,  while  CE  denotes  collisional  excitation  and  deexcitation  processes. 


X*RAY  LASER  TRANSITIONS  IN  NE-LIKE  SELENIUM 
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Figure  0.  The  4  2pB3a  and  10  2pB3p  multiplet  states  of  the  neon-like  ionization  stage  are  drawi 
to  scale  m  an  energy  diagram.  The  two  J*0-1  (denoted  A  and  C)  and  the  two  J=2  1  (denoted  B  ant 
D)  transitions  of  importance  in  neon-like  x-ray  lasers  are  identified  along  with  their  wavelengths 
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Figure  10.  The  gain  coefficients  in  CRE  that  were  calculated  using  /03  and  /*  electron 
-ustnbutions  are  shown  as  a  function  of  ion  density  for  the  J=0-1  transition  denoted  C  in  Fig.  (9). 
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Figure  11.  The  gain  coefficients  in  CRE  that  were  calculated  using  /£  and  /08  electron 
distributions  are  shown  as  a  function  of  ion  density  for  the  Je2-1  transition  denoted  D  in  Fig.  (9). 


